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Abstract 



The notion of blossom in extended Chebyshev spaces offers adequate generaUza- 
tions and extra-utihties to the tools for free-form design schemes. Unfortunately, 
. such advantages are often overshadowed by the complexity of the resulting al- 

■^ ' gorithms. In this work, we show that for the case of Miintz spaces with integer 

^^ , exponents, the notion of Chebyshev blossom leads to elegant algorithms whose 

complexities are embedded in the combinatorics of Schur functions. We express 
the blossom and the pseudo-afhnity property in Miintz spaces in term of Schur 
functions. We derive an explicit expression of the Chebyshev-Bernstein basis via 
an inductive argument on nested Miintz spaces. We also reveal a simple algo- 
rithm for the dimension elevation process. Free-form design schemes in Miintz 
spaces with Young diagrams as shape parameter will be discussed. 

> ■ 

fvj ' Keywords: Extended Chebyshev systems, Chebyshev blossom. Computer 
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t^^ ' 1. Introduction 

o. 

Representing a polynomial on an interval by its Bezier points is a common 
practice in the field of computer aided geometric design l5i] . Namely, a polyno- 
mial F of degree n, can be written as 

n 

F{t) = Y,B7{t)P..- t^[aM 

i=0 

where i?f(i) = (")a(i)'/3(i)"~*, i = 0, ...,n is the Bernstein basis in the space 
Vn of polynomials of degree n with respect to the interval [a, 6], a{t) and /3(t) 
are the barycentric coordinates of the point t with respect to the interval [a, 6], 
i.e., a{t) = {t — a)/{b — a) and (3{t) = [b — t)/{b — a). Important features of 
this representation are that the piecewise linear interpolant of the Bezier points 
{Po,Pi, ...jPn) reflects, to a certain extent, the shape of the polynomial curve 
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and that the end-segments [Po,^i] and [P„_i,P„] are tangents to the curve 
at the point F{a) and F{b) respectively. Furthermore, the curve hes in the 
convex huU of the control polygon and the de Casteljau algorithm leads to an 
efficient method for the evaluation of the polynomial from its control points. 
The total positivity of the Bernstein basis gives rise to many shape preserving 
properties. For example, the diminishing variation property ,1] ensures that the 
number of times an arbitrary hyperplane crosses the curve is no more than the 
number of times that hyperplane crosses the control polygon. The notion of 
blossom introduced by Ramshaw jl4 | offers an elegant and unifying approach 
to the understanding of the many aspects of the theory of Bezier curves. The 
fundamental idea of blossoming is that for any polynomial function F of degree 
n there exists a unique function /(ui, U2, ..., u„) that is n-affine (i.e., / is a poly- 
nomial of degree less than or equal to 1 with respect to each separate variable) , 
symmetric (i.e., /(mi,-U2, •■•,-"«) = fiua{i),u„{2),—,Ucrin)) for any permutation 
(7 on {1, 2, ..., n}) and satisfies /(i, i, ..., t) = F{t) for every i G R. The function 
/ is called the blossom or polar form of F. The control points of the polynomial 
F with respect to the interval [a, b] are then expressed in term of the blossom 
as Pi = f{a^~^,b^) for i = 0, ...,n. The multi-affinity of the blossom leads in 
a very natural way to the de Casteljau algorithm. Moreover, the blossom of a 
polynomial has a simple expression, namely if the polynomial F is expressed in 
the monomial basis as F{t) — J27=o '^^^^ then its blossom is given by 

n 
f{ui,U2,-.,Un) =^-^efc(Mi,M2, •••,-««), 

where ek{ui, U2, .-., Un) is the kth elementary symmetric function in the variables 
ui, ...,Un- The concept of blossoming was extended by Pottmann [i^ to include 
any linear space E = span (1, 01, 02, ..., 0„) such that span {(f>[, 02, ..., ^'J is an 
extended Chebyshev space of order n on an interval. The proposed extension 
bears striking similarities to the polynomial framework and in which notions of 
control points, de Casteljau algorithm and generalized Bernstein basis can be 
defined. Moreover, the emergence of the interval of interest as a shape parameter 
makes this extension fundamental in free-form curve design. However, while the 
expression of the blossom in the space of polynomials is simple, the expression of 
the blossom in a generic extended Chebyshev space is much more complicated in 
general. Thereby, leading to more complicated subdivision schemes. The main 
objective of this paper is to show that at least for the case of Miintz spaces with 
integer exponents, the notion of blossoming provides us with an elegant theory 
in which the resulting algorithms could be understood and made easy once we 
invoke the notion of Schur functions. The paper is organized as follows : In the 
second section, we review the basic properties of Chebyshev blossoming [8|, |9| . 
Emphasis will be given to the notions that will be needed within this work, such 
as the definition of Chebyshev blossom, the pseudo-affinity property, characteri- 
zations of Chebyshev-Bernstein basis and the process of dimension elevation. In 
section 3, we review the fundamentals of Chen iterated integrals. Introducing 
Chen iterated integral has a twofold aims. Firstly it will allow us to give an 
interesting determinantal expression of the Chebyshev blossom of Chebyshev 
functions defined in terms of the so-called weight functions [12], thereby, al- 
lowing the main result of the section to be used in different contexts than the 
one of Miintz spaces. Second, the determinantal expression will provide us, in 



section 5, with the Chebyshev blossom of Miintz spaces with integer exponents 
without resorting to solving linear systems. The relevant properties of Schur 
functions will be recalled in section 4. In section 5, we give the expression of the 
Chebyshev blossom of Miintz spaces with integer exponents in terms of Schur 
functions. The main result of this section is essentially the same as in [ICI, 11] in 
which a different convention was adopted and the connection with Schur func- 
tions seems to not to have been noticed. Several fundamental examples that will 
guide us throughout this work will be given. Using the Dodgson condensation 
formula, an expression of the pseudo- affinity property in terms of Schur func- 
tions will be given in Section 6. Such expression will be fundamental, through 
section?, in deriving an explicit expression of the Chebyshev-Bernstein basis in 
any Miintz space with integer exponents. Note that there is only a single case in 
which an explicit expression of Chebyshev-Bernstein basis is known, namely the 
space span{l,i'^"'"^,i'^+^,i'"'+^, ...,i'^+"} where A; is a positive integer [IjJ. Our 
strategy for deriving such an explicit expression consists of two steps. First, we 
show that, although the de Casteljau algorithm is not able to provide us with 
meaningful expressions of the Chebyshev-Bernstein basis, it will allow us to gain 
extra information on the derivatives of these bases. Then, the explicit expres- 
sion will be obtained via a dimension elevation process and some combinatorial 
manipulations on nested Miintz spaces. The Chebyshev-Bernstein bases can be 
defined without resorting to the notion of Chebyshev blossom. Therefore, this 
section shows, in particular, the importance of the notion of Chebyshev blos- 
som in solving this specific problem. In section 8, we give a simple algorithm 
for the dimension elevation process. The idea of using Young diagrams as shape 
parameter for free-from design and for the problem of continuity of composite 
Chebyshev-Bezier curves will be discussed. Expression for the derivative of the 
Chebyshev-Bernstein basis will also be given. 

2. Chebyshev blossom and Chebyshev-Bernstein basis 

In this section, we review the basic properties of Chebyshev blossoming and 
provide the relevant informations that will be used within this work. We will 
mainly follow the terminology and the notations of the excellent report [8j. 
Although, there is optimal smoothness conditions on Chebyshev functions in 
order to define the Chebyshev blossom, we will assume here, for simplicity, that 
all the functions that we encounter are infinitely differentiable. 

Chebyshev blossom: Let / denote a non-empty real interval, and let <j) — 
(01, ip2, ■■■,(f>n)'^ be a C°° function from the interval / into M" (the space M" is 
viewed as an n-dimensional affine space). Let us assume that the linear space 
D£((j)) = span(<f>\^<f>'2, ■■■,4>'n) i^ ^^ n— dimensional extended Chebyshev space 
on /, i.e., each non-zero element of this space vanishes (counting multiplicities) 
at most n — 1 times on /. In this case, we say that the function is a Chebyshev 
function of order n on I. The linear space £(</)) = span{l^(j)i,(f>2, ...,0„) is an 
(r7,-|-l)-dimensional extended Chebyshev space that we call the Chebyshev space 
associated with the Chebyshev function (/). If for any real number t in J, we 
denote by OsCi(p{t) the osculating flat of order i of the function </) at the point 
t, i.e., 

Osc,</)(t) = {</>(t) -I- ai</>'(t) + ... + «,</)« (t) I ai,...,a, gR}, 



then the assumption that </) is a Chebyshev function of order n imply that for all 
t € I and for all i = 0, ..., n, the osculating flat OsCi4>{t) is an j-affine dimensional 
space |8|. Moreover, it can be shown that for all distinct points ti, ..., r^ in the 
interval / and all positive integers /ii, ...Tfir such that X]fc=i t^k = "ni < n, we 
have 

dim n[_i OsCn-i_i^(t){Tk) ^ n-m. (1) 

In particular, if in equation ([T|) we have m — n, then the intersection consists 
of a single point in R", which we label as (^(t{'\ T2^ , ...,t^''), i.e.. 

The previous construction provides us with a function (p = ((/3i, (/?2, .-., t/^n)"^ 
from /"■ into K" with the following straightforward properties: The function (p 
is symmetric in its arguments and its restriction to the diagonal of /" is equal 
to (j) i.e., ip{t,t, ...,t) = (j){t). The function (p is called the Chebyshev blossom 
of the function 0. Note that the definition of the Chebyhev blossom imply in 
particular that if we are given n pairwise distinct real numbers Ui , i = 1 , . . . , n in 
the interval /, then the Chebyshev blossom value X = (p{ui, ...,m„) is given by 
the solution of the linear system 

det (X -(/)(«,), 0'K), ■•■,0^"~''(w»)) =0, i^l,...,n. (2) 

The pseudo-affinity property: Another fundamental property of Cheby- 
shev blossom is the notion of pseudo-affinity. Let assume given (n — 1) real 
numbers T = (ui, U2, ..., Un-i) = (tCSt^^, ..., r;^"^) (^i + ... + ^r = n — 1) in the 
interval /. According to equation ([T]), the affine space 

is an afhne line. Therefore, for any t in the interval /, the point 

ip{ui,...,Un-l,t) 

belongs to the line L. In other word, there exists a function a such that for any 
distinct numbers a and b in the interval /, and for any i e /, we have 

(p{ui,...,Un-i,t) = (1 ~ a{t)) ip{ui, ...,Un-i,a) + ait)ip{ui,...,u„^i,b). (3) 

Moreover, it is shown in [8| that the function a is a C°° strictly monotonic 
function from the interval / to R satisfying a{a) — and a(b) = 0. The function 
a will be called the pseudo-affinity factor associated with the Chebyshev space 
£{(/))■ In general, the function a depends on the interval [a, 6], the real numbers 
Ui,i = 1, ..., n — 1 as well as the parameter t. To stress this dependence, we will 
often write the pseudo-afhnity factor as a{ui, ..., u„_i; a, 6, t). 

Chebyshev-Bernstein Basis: Given two real numbers a and b in the interval 
I {a < b), and denote by 11^;, k = 0, ...,n, the (n + 1) points defined as 

n, = ^(a"-\6'). 

The points Hi are affinely independent in R" [8]. Therefore, there exist {n + 1) 
functions B^, k = 0, ...,n such that for any f e / 

n n 

m^Y.Bkim. and ^Sni) = l. 

fe=0 fe=0 



The functions Bq ,..., B"^ , ..., B" form a basis of the Chebyshev space £(</>), called 
the Chebyshev-Bernstein basis of the space £{(!)) with respect to the interval 
[a, b]. In this work, we will use the following characterization of the Chebyshev- 
Bernstein basis [9[: 

Theorem 1. The Chebyshev-Bernstein basis (_Bq , ..., S^, ..., B") with respect 
to the interval [a, b] C I, is the unique normalized basis of the space £{4>) such 
that for fc = 0, ..., n, B^ vanishes k times at a and n — k times at b. 

£^(0)-functions and its blossom: A function F from the interval / into R™, 
TO < n is called a f ((/))-function if all its components belong to the space £{(j)) 
i.e., there exists an affine map h on R" such that F = h o (j). The Chebyshev 
blossom of F is then defined as the affine image of the Chebyshev blossom of 
under the map /i, i.e., f = h o ip. We define the Chebyshev- Bezier points with 
respect to an interval [a, b] C / of a f ((/))-function F by 

P, = /(a"-\fe'), z = 0,...,n. 

As the Chebyshev blossom / of the function F inherits the pseudo-affinity prop- 
erty ([3]), the value F{t) can be computed as an affine combination of the points 
Pi,i = 0, ...,n, leading to the so-called de Casteljau algorithm. Note, also that 
the function F can be written as 



Fit) = j2Bum, t^i, 



k=0 

where (Bg , ..., i?Jf, ..., _B") is the Chebyshev-Bernstein basis of the space £{(j)) 
with respect to the interval [a, b]. 

Dimension elevation process: Consider another Chebyshev function i/j of 
order n + 1 on the same interval / and such that £(0) C £{iIj)- Let F be £{(j>)- 
function and denote by Pi, « = 0, ..., n its Chebyshev-Bezier points with respect 
to the interval [a, b]. The function F can also be viewed as an f (■(/;)-function and 
then having different Chebyshev-Bezier points Pi,i — 0, ...,n-\-l with respect to 
the interval [a, b]. From the definition of the Chebyshev blossom we necessarily 
have Pq = Pq and P„+i = P„. Moreover, it can be shown [Sj that there exist 
real numbers ^i g]0, 1[, i = 1, ..., n such that 

P^ = {l-^r)P^-l+^,P^, i = l,...,n. (4) 

3. Chen Iterated Integrals and Chebyshev Blossom 

In this section, we give an expression of the Chebyshev blossom of Chebyshev 
functions defined in term of the so-called weight functions. The notion of Chen 
iterated integrals [J] and their properties reveal to be fundamental in deriving 
such expression. Due to the simplicity of the proofs of the properties of Chen 
iterated integrals and for the sake of completeness, we will include such proofs 
in this section. 

Chen iterated integral: Let wi, W2, ••■, w„ be C°° functions on a non-empty 
real interval /. Let a and b two real numbers in /. The Chen iterated integral 
is defined iteratively as follows : 

Lt'^ - / u^iit)dt, 



and for p = 2, ..., n, we define 

J a 

Therefore, tlie Chen iterated integral can be written as 

Lt;^X...u.^^ I Wl(*l) / ^2(^2) / .... / UJnitn) dt„dt^_i...dti, 

Ja Ja Ja Ja 

or if a < & as 

4°;i...^„ = / UJl{tl)uj2{t2)...0Jn{tn)dh...dt„, (5) 






where A„ is the n-simplex in E" 

A« = {(ii,i2,...,irO eM" I 6>ti >t2 > ... >t„ >a}. 
Chen iterated integrals have the following properties [3|, Ij] 
Proposition 1. For any real numbers a, b and c in the interval I , we have 

jWM _ (_'i\nT[b,a] ic\ 

and 

n 

T[a,b] —\^ /"['='''] t\°-'C] (y\ 

with the convention that iifi^a-.w^ = 1 if r = or r > n. 

Proof. To prove ^, we can proceed as follows. We first remark that 

pb pti pt2 ptn^i 

Lt;X...u:^ = / Wi(ii) / UJ2it2) / .... / Wn(in) dt^dtn-l-dt^ 

Ja Ja Ja Ja 

h ph pb ph 

LO„{tn) I LO.n-l{tn-l) I ■■■■ I LOi{ti) dtidt2...dtn. 

Jtru •'tn-2 Jt2 

Then, we switch the limit of integration at each level starting from wi. Equation 
([7]) can be proven by induction on the number of weight functions iOi. The 
equality is obvious for n = 1. Let us assume the equality to be true for the {n — 1) 
weight functions uj2, ..., w„. Replacing 6 in ([7]) by a variable t and differentiating 
both sides of the equation with respect to t shows, by the induction hypothesis, 
that there exists a constant K such that 

n 

j^[a,t] ^\^ Li'^^t] il-^.c] +K (8) 

i=0 

Taking i = c in the last equation, shows that the constant K — 0. D 

If we take a = 6 in ([7]), and taking into account ([S]), we arrive, after renaming 
the variables, to the following 



Corollary 1. For any a and b in the interval I , we have 



YS-^"'^^-^^^r 



b] ^ T [a,b] 

...UJi^l LJnUJn — i ...UJi' 



(9) 



Remark 1. Probably the most imporant property of Chen iterated integrals is 
the so-called shuffle product of two Chen iterated integrals [3', 4] . In our present 
context of Chebyshev blossom in Miintz spaces, such property is not needed. 
However, in a future contribution we will exhibit it importance for Chebyshev 
blossom of Chebyshev functions defined in terms of weight functions. 

A determinant formulas: Let 101,102, ■■■,^n be C°° functions on a real 
interval [a,b]. Denote by /['''''l(a;i,a;2, ...,a;„) = (aij)i<ij<n the square matrix 
of order n given by 



fl" = -^i?,' 7 ai,i+i = 1. i^l,...,n, 



and 



a,j =0 if j>i + l, a,;j = iL";!,Ui--. if ■?' < *■ 
The matrix /["'^l(a;i,CL'2, ■■■,ijJn) has the form 



j-[a,b] 


1 













\ 




.[a,b] 


1 





0... 







j-[a,b] 


j-la,b] 

l^UJ2'^3 


J^U3 


1 


0... 








[aM 



[a,b] 



[a,b] 



\ T l":"I r l"."l r I' 



(10) 



T [o-M J- [a,b] I 



Proposition 2. The determinant of the matrix B'^'"{uji, 102, ■■■,U!n) is given by 
det(/[«-^lK,^2,...,^„)) = LL?i„_,...c.,- (11) 

Proof. We proceed by induction on n. For n — 2, we have 
det(/['^^''lK,a.2)) = Ltt^Ll^;'^ - ii?;3,. 

Equation ([9]) for n = 2, gives 

T[a,b] _ T[a,b]T[a,b] _ r[a,6] 

UJ2ljJl ^1 CJ2 UJiUJ2^ 

thereby, showing (fTTj) for n = 2. Let us assume (jlip to be true for all m < n. 
Now, by expanding the determinant B'^'"{u!i,uj2, ...,a;„) down the first column, 
we shall obtain 



det(/['^-^lK,^2,...,^„)) = ^(-l)*-iiL";S,...^,det(/['^'''l(w,+i,w,+2,...,u;„). 
By the inductive hypothesis, we then have 

n 

det(/[''^^lK,a.2,...,a.„))=^(-l)-iLL";3^...^^iL?i_,.....+,- 
Applying again Corollary [1] leads to the desired result. D 



Chen iterated integral and Chebyshev blossom: Let a;i,...,w„ be C°° 

functions non-vanishing on a real interval / and defined on an interval J D I. 
Let a be a fixed real number in the interval J, then it is well known [12l | that 
the function 

</.(t)=.(LLri,Lfci,...,ifci....„f (12) 

is a Chebyshev function of order n on the interval /. In the following, we show 
that the Wronskian of (j) has an interesting expression in terms of Chen interated 
integrals, more precisely, we have 

Proposition 3. For any real number t in the interval I , the Chebyshev function 
4> in \1'^) satisfies 

Proof. We first notice that 

cb'{t)=u;,{t)[l,L^^^;^Ll:;;l,...,Ll:;;l,,,^^f . 

Moreover, by a simple inductive argument, it can be shown that for 2 < fc < 
71 — 1, there exist differentiable functions pt^k such that 

fc-i 
^('='(i) = 5]p.,fe(t)0(')(t) +wi(iV2(t)...Wfe(t)*fe(i), 

2=1 

where ^fc(t) is given by 

fe-i 

^fcfi) ^ (0X^,1,^1"^*' ,Lt'^ u, ,-^1"'*' u, u, f- 

By noticing that for k ~ 2, ...n — 1, ^fc(t) is the {k + l)th column vector of the 
matrix /[°'*1 (wi, W2, •■•, w„) defined in (fTU)) . while the first and the second column 
of /[°'*l(a;i,W2, •••,'^n) are 0(i) and (j}'{t)/uji{t) respectively, we conclude that 

det(0(i),0'(i),--,0"~'W)=<"'W^r'W-^n-iWdet(/['^-*lK,^2,...,c^„)). 
The proof then results from Proposition [51 D 

Let us define the following set of functions $i, i = 1, ..., n by 

$1 (t) = fil"'*] L^"'*! l'"'*) 

and for i = 2, ..., n — 1 

^fi") = fL[°^*l L[°'*1 L^"'*! L'"'*! 

and 



T 



We have 



Proposition 4. For i = l,...,n — 1, we have 

and 

det($;(i),...,$(:-i)(t))=<-i(t).....„_i(i). 

Proof. We will show the proposition by induction on the index i. Let us start 
with the determinant formula for $i. We have <J?i(i) — aji(t)ili(i), where the 
function J7i is given by 



Therefore, we have 

Applying Proposition [3] to the Chebyshev function Vti gives 

dct(f]i(i),fii(i),...,f7i"-'^)-a.r^(i)...a.„-i(t)ii?i„_,...^,. 

Therefore, we have shown the proposition for $i. Let us assume the proposition 
to be true for any j such that 1 < j < i- We have $^(t) — uji{t)fli{t), where fti 
is given by 

n-(f) — (1 r["'*l r[«'*l ri".*] rl<^'*] )'^ 

Therefore, we have 

det($Ki),-,$!""'^W) ^ior\t)det{n,{t),n[{t),...M['"'^) ■ (13) 

Expanding the determinant det (i7i (i) , 57^ (i) , . . . , i7^" ) down the first row, shows 
that det{n^{t),n'^{t), ...,r2|""^') = dct{p[{t), ...,p^""^^(i)) where pi is given by 

01 m = fil'"'*! l["'*1 il'^'*! lI'^'*! F 

By the induction hypothesis, we have 

Inserting the result of the last equation into ([T3| leads to the desired result. D 



Let (p be the Chebyshev function of order n on an interval / defined in p^ . Let 
us denote by 0*(i) the function 

0*(i) = {^iit)^<j^m,...,K^,{t)f = {l^::'U,l^::'U,lI^:!\i] 

Using the notation 

D{fi,...,fn;xi,...,Xn) = det {fj{xk)); l<j,k<n, 
we have the following expression of the Chebyshev blossom of the function 



Theorem 2. For any pairwise distinct real numbers wi,...,w„ in the interval 
I, the Chebyshev blossom of the function (p is given by (p ~ {(pi, ...,ipn)'^ , where 
ifii is given by 

D{(j)l,...,(j)*,(p*^2,...,(p*^^j^;ui,...,Un) 
(/3,(mi,...,w„) = — ^ . (14) 

Proof. From ©, a point X — {Xi, ..., Xn)'^ in R" belongs to the intersection of 
the osculating flats of order n — 1 at the points 4>(ui) if and only if X satisfies 
the linear system 

det(X, 0'K), ..., 0*""'^ {u^) = det(0(u,), 0'K), ..., <^'""^) M) 1 = 1,2, ..., n. 

Using Proposition |3] and Proposition HI the last linear system can be rewritten 
as 

n-l 

Therefore, the statement of the theorem is nothing but the Cramer rule for 
solving linear systems. D 

If in Theorem[2]some of the real numbers Ui coincident, then we can compute 
the Chebyshev blossom from ([T^ by a straightforward iterative application of 
the I'Hopital's rule. 

4. Young Diagrams and Schur Functions 

In this section, we fix notations and review some basic concepts in the theory 
of Schur functions. We will follow the standard Macdonald's notations [7[ 

Schur functions: A sequence of non-increasing non- negative integers 

(Ai,A2,...,A„...), Ai >A2>... >A, >... (15) 

containing only finitely many non-zero terms is called a partition. The total 
number of non-zero components, ^(A), is called the length of the partition A. 
We will always ignore the difference between two partitions that differ only in 
the number of their trailing zeros. The non-zero A^ of the partition in (jlSp 
will be called the parts of A. The weight |A| of a partition A is defined as the 
sum its parts i.e., |A| = X^tli ^i- We will find it sometimes convenient to write 
a partition by the common notation that indicate the number of times each 
integer appears as a part in the partition, for example we write the partition 
A = (4,4,4,3,3, 1) as A = (4'^,3^, 1). Given a partition A, the Schur symmetric 
function 5*^(^1, ..., u„), where n > 1{X) is an element of the ring Z[wi,...,u„] 
defined as the ratio of two determinants 

c r \ det(u/"''"~-')i<i.j<„ 

det(Uj ■')i<i,j<n 

The denominator on the right-hand side of (|16p is the Vandermonde determi- 
nant, equal to the product 

V{ui,...,Un) ^ Y[ (u^-Uj). 
l<?'<j<n 

10 



We will adopt the convention that Sx{ui, ...,Un) = if /(A) > n. From the 
definition, the Schur function associated with the empty partition A = (0, ..., 0, ..) 
is S\{ui, ...,Un) = 1. For the partition A = (r), the Schur function S\ is the 
complete symmetric function hr i.e., 

S(r)iui,U2,...,Un) ^ hr{ui,...,Un) ^ ^ Ui^Ui^.-.Ui^, 

while for the partition A = (f) with r < n, the Schur function S'(ir) is given by 
the elementary symmetric function e^ i.e, 

S(ir){ui,U2,...,Un) ^ er{ui,...,Un) ^ ^ Ui^Ui^...Ui^ . 

il<i2<---<ir 

A direct consequence of the definition is the following 

'5(Ai + ia2 + l,...a„ + l)(Ml, ..., U„) = MlM2-.-MnS'(Ai,A2,...a„)(Ml, ..., U„). (17) 

The Schur function 5a can be expressed in terms of the complete symmetric 
functions through the Jacobi-Trudi formula 



S\ = det(/iAi-i+j)i 



^i.j<:n ' 



(18) 



where we assume that hm = if to < 0. The conjugate, A', of a partition A 
is the partition whose Young diagram is the transpose of the Young diagram 
of A, equivalently A^ = Card{j\Xj > i}. Using the conjugate partition, the 
Schur function can be expressed in term of the elementary symmetric functions 
through the Nagelsbach-Kostka formula 



5a = dct(eA;-i+j)^ 



<i,j<n ' 



(19) 



where we assume that e^ = if m < 0. Throughout this work, we will use the 
notation 

to mean the evaluation of the Schur function in which the argument ui is re- 
peated TO-i times, the argument U2 is repeated TO2 times and so on. 

Combinatorial definition of Schur functions: The Young diagram of a 
partition A is a sequence of 1{X) left -justified row of boxes, with the number of 
boxes in the ith row being A^ for each i. A box x = {i,j) in the diagram of A 
is the box in row i from the top and column j from the left. For example the 
Young diagram of the partition A — (5, 4, 2) and the coordinate of its boxes are 



A =(5,4,2) 



(1,1) 



(2,i; 



(3,i; 



(1,2) 



(2,2) 



(3,2) 



(1,3) 



(2,3) 



(1,4) 



(2,4) 



(1,5) 



A semi-standard tableau T^ with entries less or equal to n is a filling-in the 
boxes of A with numbers from {l,2,...,n} making the rows increasing when 
read from left to right and the column strictly increasing when read from the 
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top to bottom. We say that the shape of T'*' is A. For each semi-standard 
tableau T^ of the shape A, we denote by pi the number of occurrence of the 
number i in the semi-standard tableau T^ . The weight of T^ is then defined as 
the monomial 

For a given partition A of length at most n, the Schur function Sx{ui, ..., w„) is 
given by 



Sx{ui,U2,...,Un) = ^'i 



where the sum run over all the semi-standard tableaux of shape A and entries 
at most n. 

Example 1. Consider the partition A — (2, 1) and n = 3. Then, the Young 
diagram of A and the complete list of semi-standard tableaux of shape A are 



1 


1 


1 


1 


1 


2 


1 


3 


1 


2 


1 


3 


2 


2 


2 


3 


2_ 




_3_ 




1 




2_ 




2_ 




_3_ 




J_ 




i_ 





Therefore, the Schur function associated with the partition A is given by 

Sx{ui,U2,U3) = u\u2 -f M1M3 + 2uiU2'U3 -|- U^U^ + U2U\ + Uiu\ + Uiu\. 

Giambelli formula: The Young diagram of a partition A is said to be a hook 
diagram if the partition A is of the shape A = (p -I- 1, f) i.e., 

P+i 



In Frobenius notation, we write the partition A as {p\q). Expanding the 
Jacobi-Trudi formula (J18p along the top row, shows that the Schur function 
associated with the partition {p\q) is given by 



S{p\q) — hp+lCq — /lp+2eg-l + .... + (— l)'/lp+q+l. 

Any partition A can be represented in Frobenius notation as 



A 



[ai, 



, ar|/3i,...,/3r). 



(20) 



(21) 



where r is the number of boxes in the main diagonal of the Young diagram of 
A and for i = 1, ...,r, a^ (resp. /3i) is the number of boxes in the ith row (resp. 
the «th column) of A to the right of {i,i) (resp. below («,«)). For example the 
partition A = (6,4, 2, 1^), depicted below, can be written in Frobenius notation 
as A= (5,2|4, 1) 
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With the decomposition (pij) of A in hook diagrams, the Giambelh formula 
states that 

Sx = dct{S(^ai\l3j))l<i,j<r (22) 

We win adopt the convention that S(^a\i3) = if a or /3 are negatives. 

Hook length formula: The hook-length of a partition A at a box x — {i,j) 
is defined to be h{x) — Xi + X'^ — i ~ j + 1, where A' is the conjugate partition of 
A. In other word the hook-length at the box x is the number of boxes that are 
in the same row to the right of it plus those boxes in the same column below it, 
plus one (for the box itself). The content of the partition A at the box x — {i,j) 
is defined as c{x) = j — i. The hook-length and the content of every box of the 
partition A = (5, 4, 2) is given as 



MA) = 



7 


6 


4 


3 


1 


5 


4 


2 


1 




2 


1 









Content{X) — 






1 


2 


3 


4 


-1 





1 


2 




-2 


-1 









With these notations, the number of semi-standard tableaux of shape A with 
entries at most n is given by the so-called hook-length formula as 



h{n) - 5^(1X71) = n 



n + c(x) 
h{x) ' 



In particular, we have the following useful hook-length formulas 

1' 



and 



/(!-)(") = 

/(Pk)(") 



f{r){n) 



n fn + p 
P+q+l\ P , 



+ r- 
r 

n — 1 

q 



(23) 



(24) 



(25) 



We will adopt the convention that for every integer n, the hook-length of the 
empty partition A = (0, 0, ...) is given by /0(n) = 1. 

Skew Schur functions and Branching rule: Given two partitions, A and 
/x, such that // C A i.e., m < Xi, i > I, a, Young diagram with skew shape X/^ is 
the Young diagram of A with the Young diagram of /i removed from its upper 
left-hand corner. Note that the standard shape A is just the skew shape X/fi 
with /i = 0. For example, we have 



(4,3,1)7(2,1) 



The skew Schur function S\/^ is defined as 

Sx/t,{ui,U2,-.,Un) ^ ^ X 



rpX/ fj. 



where the sum run over all the semi-standard tableaux of shape X/fi and entries 
at most n. Skew Schur functions have a determinant expression as 

^x/p, — det{hxi-i_ij-i+j)i<ij<n 
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Using the skew Schur functions, we have the foUowing branching rule 

Sx{ui,...,Uj,Uj+i,...,Un) = ^ S^{ui, ...,Uj)S),/^{Uj+i, ...,Un). (26) 

Particularly interesting for this work, the following two branching rules 

5a(ui,...,u„-i,w„) = ^ 5^(wi,...,m„_i)mIj^I"I^I, (27) 

where the sum is over are the interlacing partitions fi i.e., partition fj, — (/ii, ..., 
fj,n~i) such that 

Ai > Ml > A2 > ...Ai„-i > A„, (28) 

and 



Ai 



S\{ui, ...,u„_i,u„) = ^S'a/(j)(ui,...,u„_i)m^. 



(29) 



5. Blossom in Miintz space with positive integer pow^ers 

It is well known that for any positive real numbers < si < S2 < ... < Sn, 
the function (f>{t) — {t'^'^^t^^, ...,i''")"^ is a Chebyshev function of order n on the 
interval ]0,cxd[. In this section, we give the Chebyshev blossom of the function 
(j) in case the parameters Si,i — l,...,n are positive integers. We will first 
associate the sequence (si, ...,s„) with a partition A that will allow us to give 
the expression of the blossom in terms of Schur functions. We will first start 
with a definition 

Definition 1. Let A = (Ai, A2, ..., A„) be a partition of length /(A) at most n. 
The Miintz tableau associated to the partition A is given by a sequence of n+ 1 
partitions (A(°\ A^^^ A^^^ ..., A^"^) defined as follows: 

AW = (A2,A3,...,A„), 

for j = 1, 2, ...n — 1 

AW = (Ai + 1, A2 + 1, ..., A, + 1, A,+2, ..., A„) 

and 

A(") = (Ai + 1,A2 + 1,...,A„ + 1). 

To remember the construction of the Miintz tableau we can remark that the 
partition A*-"' is obtained form the partition A by deleting the first row. The 
partition A*-*^) will play an important role in this work and will be called the 
bottom partition of X. The partition A*-*-* is obtained by adding a box to the first 
i rows of the partition A, deleting the i + 1 row and keeping all the other rows 
the same. 

Example 2. The Miintz tableau associated with the partition A = (4, 2) and 
n — 'i is depicted as 



A 



A(o) 



Ad 



A(2 



A(3) = 
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To a given partition A — (Ai,...,A„) of length at most n, we define the 
following Chebyshev function of order n 

(t){t) = ft^i-^2 + lj-Xl-\,+2^^_^^\i-\„ + {n-l)^^\l+n\ (3Q) 

The associated Chebyshev space £{(j)) will be denoted by Sxin) and will be 
called the Miintz space associated with the partition A. The function will be 
called the Miintz function associated with A and conversely, the partition A will 
be called the partition associated with the function (f>. We have the following 

Theorem 3. For any sequence (ui, U2, ■■•, Wn) g]0,+(X)[", the blossom ip — 
{(pi,(p2, ■■■,^n)'^ of the Chebyshev curve (f> given in i30\) is given by 

I N f\(o){n)Sx(i){ui,U2,...,u„) 
ip,{ui,U2,...,U„) = -— -, 

/A(i) [njoxm (ui,U2, ..., Un) 

where (A^*^-*, A*-^', ..., A*-"') is the Miintz tableau associated with the partition A 
and fp,{n) refers to the number of semi- standard tableaux of shape fi and entries 
at most n. 

Proof. We first assume that all the positive real numbers Ui,i — 1, ...,n are pair- 
wise distinct. Consider the functions aji,aj2, •••, w„ such that for i = 1, 2, ...,n — 1 

0.W-4"4..... =i'^-''^^+' and 0„(t) = iL°4....„=i'^+"- (31) 

Applying successive derivatives to pip shows that there exist positive constants 
Ki,i = 1 , . . . , n such that 



W: 



(t) ^ Kit^'~^'+\ for i = l,...,n-l and ujn{t) = K^t^ 



Computing the Chen iterated integrals of the obtained function w^, i — 1, ..., n, 
shows that there exist constants Ck,k — 1, ...,n such that 

r[0,t] _ ^ j-Afc + (n-fc+l) 

From Theorem [21 the Chebyshev blossom of the function (p can be expressed as 

I \ r^' 0,et[U^ )l<i,j<n ,„„•, 

^,{ui,...,u„)=C \' J -''- , (32) 

det[Ui' )i<i,j<n 

where 

Cj', i — 1, ..., n are constants and p = (Ai + 1, ..., A^ + 1, Ai+2, ■•■, A„) = A^*-*. 

Dividing both the numerator and the denominator of the right hand side of (j32p 
by the Vandermonde determinant Y[i<i<j<ni'^i ^ ^j) leads to 

(pi(Ui,U2,...,Un)=C^- r. (33) 



Now, as the expression (|33| still make sense even if some of the Ui coincident, 
and since the process of intersecting osculating flat is a smooth process, the 
Chebyshev blossom of the function </> evaluated at any positive real numbers 
Ui,i — 1, ...,n is still given by the expression ([55]) . The value of the constants 
C,' in p3l) can be obtained as follows : From the definition of the Chebyshev 
blossom, we have (p{t,t, ...,t) — (t){t), then in particular we have ip{l, 1, ..., 1) — 
(1,1, ...,1), which gives the value of the constants as claimed by the Theorem. 

D 
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Examples section: Horizontal, vertical and hook Young diagrams occupy an 
important place in the combinatorics of Schur functions. Therefore, it is only 
natural to define the Miintz spaces associated with these particular Young di- 
agrams and carry throughout this work their fundamental properties. Some 
times, we will also give low order Miintz spaces to exhibit the use of the combi- 
natorics of Schur functions in solving particular problems. We will also define 
the staircase Miintz space as they have the particularity of being, in a sense to 
be precised, a "reparametrization" of the polynomial spaces. 

Polynomial Miintz space: Consider the Chebyshev curve of order n over 
the real line R. 

0(t) = (i,i^...,i"^. (34) 

The associated partition A is the empty partition and the space £0 (n) is the lin- 
ear space of polynomials of degree n. The bottom partition X^^' is also an empty 
partition, while the rest of the Miintz tableau is given by A^'"'^ = (l'"): fc = 1, .-., n. 
Therefore, the Chebyshev blossom oi ip — {ipi, ip2, ■■■, 'Pn)^ of the function is 
given by 

S'(ifc)(ui,...,u„) efc(Mi,...,M„) 

Combinatorial Miintz space: Consider the Chebyshev function (f){t) — 
[t.t^.t'^Y' of order 3 over the interval ]0,cx)[. The partition A associated with 
the curve is given by A = (1,1,1). The Miintz tableau associated with A is 
given by (A^ = (1, 1), A^) = (2, 1, 0), A(2) = (2, 2, 0), A(3) = (2, 2, 2)) 



A = 



— 


A(o) = 


— 


aW = 




A(2)- 




— 


A(3) = 







Therefore, the blossom ip = [fi^ 'fi-, "y^s)"^ of the function </> is given by 

/a(0)(3) 5'A(fc)(Ml,M2,M3) 



V3fc(ui,U2,M3) 



/A(fc)(3) S'a(0)(ui,U2,U3)' 



We can now proceed by computing the Schur functions associated with the 
partitions in the Miintz tableau. For the partition A*^"-* = (1, 1), we have 

S';^(0)(ui,U2,-U3) = 62(^1, W2, -"3). 

The Schur function associated with the partition A*^^^ — (2, 1) has been already 
computed in Example 1. The semi-standard tableaux associated with the par- 
tition A^^-* = (2, 2) and entries at most 3 are given by 



1 


1 


2 


2 



1 


1 


2 


3 



1 


1 


3 


3 



1 


2 


2 


3 



1 


2 


3 


3 



2_^ 
3 3 



Therefore, we have 

'5'(2,2)('"lj'«2,'"3) = U1W2 + ulu2U3 + ulul + U^UiUs + U3M2W1 + ulul- 

For the partition A*^'^-' = (2, 2, 2), we can use ([T7| to deduce that 



2„,2„,2 



S'(23)(wi,M2,W3) = UiU2U3S(i3){ui,U2,U3) = UiU2U3e3{ui, U2, U3) = UiMgU 
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Therefore, the blossom ip of the function (j) is given by 



ip{ui,U2,U:i) 



1 



[U1U2 + M1U3 + U2W3) 



3("1 +"2)("1 +"3)("2 +"3) 
4 ((,,2„| + „2„2 _,_ „2„2) + „^„2„3(„j + ^^ + ^3)) 



Elementary Miintz spaces: Let k and n be two positive integers such that 
1 < k < n. Consider the Chebyshev curve of order n over the interval ]0, oof 
defined for fc 7^ 1 by 

0(i) = (t,i^...,i^-^t'=+^...,^+lf, 

and (f){t) ~ (t^, t'^, ...,i"+^)-^ for k = 1. The partition A associated with the 
function is given by a vertical Young diagram with k boxes, i.e, A — (1*^). 
For this reason, we will call the curve <f> the kth. elementary Miintz curve and 
the space EiikAn) the fcth elementary Miintz space. The bottom partition 
^(0) ig given by A*^°-' — (l*^"^), with an associated Schur function given by 
S\{o){ui, ...jUn) — efc_i(Mi, ...,M„). The other partitions in the Miintz tableaux 
are given by 

^W ^ (•2M'=-^-i) for i = l,...,fc-l and A^'^ = (2^ l*"*^) for i^k,...,r 

The Young diagram of the partitions in the Miintz tableau are of the form 



A(o) = (fc - 1) < 



For i = l,...,fc — 1, we have A*^*^ == i 



{k-i-l) 



and for i = k, ..., n A^*) = k 
{i-k){ 



The conjugate of the partition A*-'-' for i — 1, ..., n are then given by 

AW' = (fc-l,i) for i = l,...,k-l and AW' = (i,fc) for i = k,...,n. 

Therefore, Using Theorem|3]and the Nagelsbach-Kostka formula (|19p . the Cheby- 
shev blossom ip = {ipi, (p2, ■■■,'Pn)'^ of the function is given by 



fi 



\k~l/ 



Bk-lSi - ei-iBk 



Cl)il)-L-I)il) 



efe-1 



for i = l,...,fc — 1, 



and 



y^i 



Kk-lJ 



Gke-i - ei+iCk-i 



a)(:)-G:j(.^) 



for i = k, ..., n. 
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Of a particularly interesting form is the last component of ^ as we have 

(p„(ui,...,u„) 



(fc-i) (nr^i'»i)efc(Mi,.-.,M«) 

(fe) efe_i(ui,...,u„) 



Complete Miintz spaces: Let fc be a non-negative integer and denote by (f) 
the Chebyshev curve of order n over the interval ]0, oo[ 

The partition associated with the curve (j) is given by a horizontal Young diagram 
with k boxes, i.e., A = (A;). We will call the function the fcth complete Miintz 
function and the associated space £(fe)(n) the /cth complete Miintz space. The 
bottom partition A'"-* is an empty partition, while the other partitions in the 
Miintz tableau are given by A*-*' — {k\i — 1). Therefore, the Chebyshev blossom 
ip= ((^1, ..., ipn)'^ of (j) is given by 

S(k\i-l){ui,U2,--.,Un) 

ip,{Ui,U2,...,Un) = -^ , 

/(fc|i-l)N 

where Sf^i_i\ can be expressed in term of the complete and elementary symmet- 
ric functions according to (|20l) and the normalization constants can be computed 
using equation (05|) . Note that ipi — hk+i / {^^t^) , while fn ~ e„/ife/(" ^^^ ). 

Hook Miintz spaces: Let I and n be two positive integers and let fc be a 
positive integer such that k < n. Consider the Chebyshev curve (f) of order n 
over the interval ]0, oof given by 

</)(i)-(i'+\t'+',...,i'+', *'+"+',..., i'+"+'r. 

The partition A associated with the curve (j) is given by a (Z,fc)-hook Young 
diagram, i.e., A = {l\k). Therefore, the function cj) will be called a (Z,fc)-hook 
Miintz function, while the associated space £(i\k) (n) will be called the {I, fc)-hook 
Miintz space. The bottom partition A*^°' is given by A*^"^ = (l''): while the other 
partitions in the Miintz tableau are given by 

X^'^ = il + 2,2'-\l^-') for i = l,...,fc 

and 

A(') = (^ + 2,2^1^-'=-l) for i = k + l,...,n. 

Every partition in the Miintz tableau has at most two boxes in its main diago- 
nal, thereby, making Giambelli formula (1^^ useful for the computation of the 
associated Schur functions. In Frobenius notation, the partitions in the Miintz 
tableau are given by 

X^'^ = {l + l,0\k-l,i-2) for i = l,...,k 

and 

A(*) = (? + l,0|i-l,fc- 1) for i = k + l,...,n. 

Therefore, the Chebyshev blossom ip = {ipi, ...,ipn)^ of is given by 

'P^ = -r'^^ -^-^ — -^ — - for z = l,...,fc 

/a(') ("-) efc 
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and 

where the normahzing factors fx{i) (n) can be computed using equations (j24p 
and (j25l) . In particular, we have 

(fc) 'S'(;+i|fc_i) (fc) e„5'(;|fc) 

/(/+i|/c-i)("-) efe J{i\k)\n) Ck 

Staircase Miintz spaces: Let (p — (0i, 02, . ..</>„) be a Chebyshev function 
of order n on an non-empty interval / and denote by ip its Chebyshev blossom. 
Let 9 : I — > / be a C°° strictly monotonic function. Then, the function 

= 0o(0,0,...,0) (35) 

is a Chebyshev function of order n on the interval /. Moreover, as the process 
of intersecting osculating flats is a geometrical concept depending only on the 
curve itself and not on the chosen parametrization, the Chebyshev blossom (p 
of the function (f> is given by 

^ = ipo{0,0,....,0). (36) 

For similar reasons, the pseudo-affinity factor a of the space £{(f>) is related to 
the pseudo-afflnity factor a of the space £{4)) by 

a(ui, ...,u„_i;a,6,t) = a(6'(ui), ...,d{un-i);9{a),9{b),9{t)). (37) 

Finally, if we denote by BJ^ and BJ!:, k = 0, ...,n the Chebyshev-Bernstein basis 
of the spaces £{(j)) and £^((^) respectively, then we have 

BUt) = BUO{t)), k^O,...,n. (38) 

Now, we will deal with the simplest case of a situation such ([55)) . namely, a 
reparametrization of the space of polynomials. Let I be a non-negative integer 
and consider the Chebyshev curve (f) of order n over the interval ]0, cx)[ given by 

The partition A associated with the function (j) is given by the so-called Z-staircase 
partition 

X^{nl,{n-l)l,{n^2)l,...,l). (39) 

The function (j> will be called a /-staircase Miintz function, while the associated 
Chebyshev space will be called a /-staircase Miintz space. The function cj) can 
be rewritten as 

(/.(t) - ((t(' + l))l, (t(' + l))2, ..., (iC+l))^ ..., (tC + l))")^. 

Therefore, the function is a reparametrization of the Miintz polynomial func- 
tion ([M]) . Taking the Chebyshev blossom of cj) using Theorem [3] in one hand 
and equation ([55]) in the another hand, in which d(t) = t'+^ in ([55]) . lead to a 
set of power plethysms 

ek{u[^^,u'-^^,...,u''+^) ^ /a(o) Sxm{ui,U2,-.,Un) , . 

(2) AW 5a,o,(ui,U2,...,U„)' ^ "^ 
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where (A^°\ ..., A*^")) is the Miintz tableau associated with the partition A in 
([M)) . Note that (|i(Il) is not a genuine power plethysm as we do not expand 
the quantity in the left hand of PO)) in the Schur basis. In this work, our 
interest in the staircase Miintz spaces is motivated by two facts. The first, is 
that as their pseudo-affinity factors as well as their Chebyshev-Bernstein bases 
are well known, they will play a role of reconfirming our theoretical results. 
The second fact is that, in practice, these spaces will play a sort of short-cut 
in finding explicit expressions of the Chebyshev-Bernstein basis for a generic 
Miintz space. A property of ^-staircase Young diagram that will be needed later 
is the following expression of their associate Schur functions, namely for the 
partition A given in p9p. we have 



u'+i - u[+' 



Sx{ui,U2,.:,Un,Un+l)= J^ hi{ui,Uj)= ]^ — —-^ . (41) 



Ui — Ui 
l<i<j<7i+l l<i<j<n+l ■' 



Remark 2. The definitions of elementary, complete and hook Miintz spaces in 
our previous examples depend primarily on the convention that we have adopted 
in associating a Miintz space to a partition in (PH)) . However, as it will be clear, 
once we give the expressions of the pseudo-affinity factors and the Chebyshev- 
Bernstein bases of these spaces, that the adopted convention is the most natural 
one. 

Remark 3. Theorem [3] it true even if A = (Ai, A2, ..., A„) is such that Ai > 
A2 > ..• > A„ > and A^ are real numbers. In this case, the Schur function 
should be defined only as the ratio of determinants as in (I16|) and in which we 
make use of the rHopital's rule when some or all of the arguments coincident. 
In the case the Ai are positive rational numbers, we can, in principle, write the 
associated Chebyshev function as a composition of the form psp and in which 
the Chebyshev function is associated with a true partition. For example, the 
Chebyshev function 4){t) — it's , is ^ is )^ on the interval ]0, oo[ can be written as 
0(i) = (ie , (^6)3(^6)4)^. Therefore, we can use the remarks in examples section 
related to the staircase Miintz spaces to compute the blossom of the function 0. 

Remark 4. In the proof of Theorem|21 we have decided to not to keep track of 
the exact value of the constants that naturally appears within the proof. The 
main reason for this decision is the fact that we can always use the diagonal 
coincidence property of the Chebyshev blossom to compute the final normalizing 
factors. However, if we had kept track of the constants, we would have proven 
a formula for the ratio of the hook-lengths. The fact that complete Miintz 
spaces have polynomials blossom would have then allow us to find a new proof 
for the hook-length formula (|23p for the hook Young diagrams and then using 
the Giambelli formula, we would have proven a determinant expression for the 
hook-length formula. 

In the following, we would like to draw attention that the Miintz tableau 
associated with a partition A appears naturally in the expansion of the Jacobi- 
Trudi determinant P^ . More precisely, let {}S^\ ..., A*^"^) be the Miintz tableau 
associated with a partition A of length at most n. Computing the Schur function 
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Sx(n) using the Jacobi-trudi determinant P^ by expanding the determinant up 
the last column [7], we find 

n 

'5'a(") = (^1)" 'S';,(0)^Ai+n + 2^(-l)" '/lAi + (n+l-i)5'A(i-i)- (42) 

i=2 

Dividing (1421) by Sxm and normalizing using the hook length factors f\ii){n), 
we arrive at 

Proposition 5. Let A = (Ai,...,A„) be a partition of length at most n. Let 
If — {ifi, ...jifin)'^ be the blossom of the Chebyshev function associated with the 
partition A. Then we have 

~{ ./ A(0) ("J 

where (A*^"', ..., A*^"') is the Milntz tableau associated with the partition A and 
An+l = 0. 



6. The pseudo-affinity factor 

For a given partition A of length at most n, we give an expression of the 
pseudo-affinity factor associated with the Miintz space £\ (n) in terms of Schur 
functions. The following so-called Dodgson condensation formula 6] will be 
crucial to this end. 

Proposition 6. Let A be an (n, n) matrix. Denote the submatrix of A in which 
rows ii, 12, ..., ifc and columns Ji,J2t ■■■, jk o,i"e omitted by Aj^'j^'"''f'' . Then we 
have 

det{A)det{Alf;) = dei(A})dei(A;^) - detiA'DdetiAl). (43) 

From the last proposition, we can prove the following 

Proposition 7. Let A = (Ai, A2, ..., A„) be a partition of length at most n. 
Then, for any sequence of real numbers U = (mi, W2, .-., "Un-i) cind real numbers 
X, y, we have 

{x - y)Sx{U, X, y)S^m (U) = xS\{U, x)Sx{o) {U, y) - ySx{U, y)Sx(o) {U, x), 

where X^^^ is the bottom partition of X. 

Proof. Without loss of generality, we can assume that all variables Ui,i = 
1, ...,n — l;a; and y are pairewise distinct. Let us denote by Vu the Vander- 
monde factor 

Vu = n ("«-^j)- 

l<i<j<n-l 

Now, let us apply Proposition [6] to the (n + 1, n + 1) matrix A defined as 

A = (a,,)i<,,,<„+i = x^+("+'^-^' (44) 
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where xi — x, Xi — Ui-i for i — 2, ...,n,Xn+i — y and A„+i = 0. The fohowing 
determinant formulas can be readily checked 



det{A) = {x-y)VuSx{U,x,y)\{{x-u{){u^-y)- det(A};:^) = VuS^{U) \{ u, 

n—1 n—1 

det{A\) = VuS^{U,y) J|K - y); dei(A;0 = xVuSx{U,x) J| u,{x - u,) 

i=l i=l 

n—1 n—1 

detiA'l) = yVuSx{U,y) J] u,{u, - y); det{Al) = VuS^{U,x)Y[{x - u^) 

4=1 1=1 

Upon applying (|43)) . the claim of the proposition follows. D 

At this point, we can give a Schur function representation of the pseudo- 
affinity factor as follows 

Theorem 4. The pseudo-affinity factor of the Chebyshev space £x(n) associated 
with a partition A of length at most n is given by 

t-a SxiU,a,t)SxmiU,b) 

a{U]a,b,t) = j-j-TTT — rr^ 777-7, 

b-a Sx[U,a,b)Sxm[U,t) 

and 

RfjT u +\ 1 iTT h +^ ^ b - 1 Sx{U,b,t)Sxm{U,a) 
maAt) = l-a{U^aAt) = (t—J sx{U,a,b)Sx^oAU,ty 

where U is a sequence of positive real numbers U — (ui, ..., w„_i) and X'^' is 
the bottom partition of X. 

Proof. Let (j) the Miintz function associated with the partition A, and (p = 
{ifii, ...,(pn)'^ its Chebyshev blossom. As the pseudo- affinity factor is indepen- 
dent of which £x{n) function we choose, we can work with the last component 
LPn of the blossom 

(^„(mi,...,m„) = J ^--t — . (45) 

By equation ([3]) and if we denote U = (ui, U2, ..., u„_i), we have 

trr u ^\ Vn{U,t) -ipniU,a) 

a{U;a,b,t) = — — — 777^- 

iPn[U,b) -ipn[U,a) 

Inserting pS)) into the last equation, leads to 

,.r , . ^ tSxjU, t)Sxm jU, a) - aSxjU, a)S xjq, {U, t) Sxm jU, h) 
""^ '"' ' ^ bSx{U,b)SxioAU.a)-aSx{U,a)Sx,o,{U,b)SxioAU,ty ^ ' 

Applying Proposition [7] with x — a and y ^ t to ^^ leads to the desired 
expression for the pseudo- affinity factor. Similar treatment with 1 — a leads to 
the second equation of the proposition. D 
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The pseudo-afRnity factor of the Miintz spaces defined in the examples sec- 
tion can be derived for the last proposition. For the Miintz polynomial space 
£(/t{n), the partition A and it bottom partition A^"-* are empty and therefore by 
Theorem 01 the pseudo-affinity factor is given by 

t — a 

a{ui,...,Un-i;a,b,t) = . 

b — a 

Similarly, the pseudo-afRnity factor of the kth elementary Miintz space £iik\{n) 
is given by 

/ u -,\ t- aekiui,...,Un-i,a,t)ek-i{ui,...,Un_i,b) 

a{ui,...,Un-i;a,b,t) = — . 

b- aek(ui,...,Un-i,a,b)ek~i[ui,...,Un-i,t) 

For the fcth complete Miintz space £^(j,)(r7,), we have 

t - a hk{ui,...,Un^i,a,t) 

a{ui,...,Un-i;a,b,t) = -— — . 

6- ahk[ui,...,Un-i,a,b) 

For the (Z,fc)-hook Miintz space f(/|fc)(?i), we have 

t- a S(^i\k){ui,...,Un-i,a,t)ek{ui,...,Un^i,b) 



a{ui,...,Un-i;a,b,t) 



b- a S(^i\k){ui, ...,u„„i,a, &)efe(wi,...,w„_i,t) 



Consider now the pseudo-affinity factor of the ^staircase Miintz space associated 
with the partition A in (j39p . Using the fact that 



S'a(mi,...,u„-i,w„,w„+i) =]^/l;(ui,Mj), 



Sx(<»{ui, ...,u„_i,M„) =Y\^ hi{ui,Uj), 

l<i<j<n 

and carrying out all the simplifications that appear in the computation of the 
pseudo- affinity factor, we find 

, , , t-ahi{a,t) f'+i-a'+i 

which is what is expected from the relation (1371) . 

For later use, we will need the equivalent of Proposition[71 for every partition 
A'-''-' in the Miintz tableau of the partition A. 

Proposition 8. Let \ be a partition of length at mostn and let (A'^*'\ A'^^^ ..., A'-"-') 
its Miintz tableau. Then, for any real numbers U = (mi, ...,m„_i), real numbers 
X and y, and k = 1, ...,n — 1, we have 

Sx(.o){U,x)Sxw{U,y) - Sx(.o){U,y)Sx(k){U,x) ^ {y - x)Sjj{U)Sx{U,x,y), 

where 77 is the bottom partition of X^''' i.e., 77 is the partition rj ~ (A2 + 1, A3 + 
1, ..., Afc + 1, Afc+2, •••, A„). 
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Proof. We can, without loss of generality, assume that all the variables Ui,i — 
1, ...,n — 1, X and y are pairewise distinct. Consider the (n + 1, n + 1) matrix A 
defined in pi)) . Now, construct a matrix B^ by putting the first column of the 
matrix A as the last column and putting the (fc + l)th column of the matrix A 
as the first column. The proof of the proposition is then derived by applying 
the condenstation formula P5|) to the matrix B. D 

Remark 5. Note that Proposition [5] can be used to reconfirm the fact that the 
pseudo- affinity factor associated with a Miintz space £\(n) can be computed 
from ([3]) using any component of the Chebyshev blossom. To show this fact, 
we can choose to work with the component ipk of the Chebyshev blossom and 
give the expression of the pseudo-affinity factor in a similar fashion as in the 
proof of Theorem [3] and in which this time we use the last Proposition instead 
of Proposition [T] 



7. The Chebyshev-Bernstein Basis 

The main objective of this section is to give an explicit expression in terms of 
Schur functions of the Chebyshev-Bernstein basis of the space £\ (n) associated 
with a partition A of length at most n. As the proof involve several technical 
steps, we will first give the main result and some of it consequences. We will 
explain the methodology of the proof along the coming subsections. 

Theorem 5. The Chehychev- Bernstein basis {Bq -^, B" .^, ...,i?"_,^) of the Miintz 
space associated with a partition A = (Ai, A2, ..., A„) of length at most n over an 
interval [a, b] is given by 

where _B^ is the classical Bernstein basis of the polynomial space over the interval 
[a,b] and A^*^^ is the bottom partition of X. 



To exhibit the fact that the Chebyshev-Bernstein basis B^ ^ in (1481) is indeed 
a polynomial function in t, we could use the Branching rule (j27l) as 

the sum is over are the interlacing partitions 77 i.e., partition 77 = (771, ...,?7„_i) 
such that 

Ai > 771 > A2 > ...?7«-i > A„. (49) 

Therefore, 

t^^Sx ('a"-^6^-) =Y,S,{a^^-'',b''){ab)^^^-MtM-\^"'>\, (50) 

where A'-'^' is the bottom partition of A. Any partition 7/ that satisfies P^ also 
satisfies |ry| — |A''°'| > 0. Therefore, for any fc = 0, .., n, the Chebyshev-Bernstein 
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function BJ^ x ^^ ^ polynomial function in t. We can also use the branching rule 
(HH) as 



^ ^ .7=0 



Al 



(51) 



0=0 



The term fx{n + ^)/ f\(a){n) in (j48p can be computed using the hook-length 
formula (P5)) . However, since we have a ratio of hook lengths of two related 
partitions, several simplifications will appear. In fact as the following lemma 
shows, to compute this term, we need only to form the hook length and the 
content of the first row of the partition A. 

Lemma 1. Let A = (Ai, A2, ..., A„) be a non-empty partition 0/ length at most 
n and let X-^' be its bottom, partition. Then, we have 

fx{n+l) _-pV(n + l)+CA(l,.7) 



fxm ("•) 



j=i 



hx{l,]) 



Proof. If the partition A consist of a single part A = (Ai, 0,0, ..., 0), then the 
partition A'"-' is empty and the lemma is the statement of the hook length 
formulas (1^51) . Let us assume that A = (Ai, A2, ..., As, 0, ..., 0) consists of more 
than a single part, i.e., Ai > A2 > 1. The definition of a bottom partition imply 
that for every non-empy box in the partition A'-'^', and i 7^ 1, we have 

hx{i,j) = h^m{i-l,j) and ca (i, j) = C;^(o) (? - l,i) - 1. 

Therefore, for any z 7^ 1, we have 

{n+l) + cx{i,j) n + Cx{o}{i-l,j) 



hx{i,j) 



hxm{i- l,j) 



Thereby, we have 



fx{n + l) _^{n+l) + cx{l,j) 



fxm (") 



j=i 



hx{l,j) 



n 



Example 3. Consider the partition A = (4,2,1). Then the content and the 
hook length of the first row are given by 



Content{X) = 






1 


2 3 















MA) = 



6 


4 


2 1 















Therefore, 



fxjn + l) _ {n + 1) (n + 2) {n + 3) (n + 4) 
fxmin) "6 4 2 1 ■ 
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Using Theorem [5J we can give the explicit expression of the Chebyshev- 
Bernstein basis associated with the Miintz spaces defined in the examples section 

Combinatorial Miintz space: Let {B^ ,^ ^-.^ Bf ,^ 2)1 (^2 (2 2)' -^3 (2 2)) ^^ ^^^ 
Chebyshcv-Bcrnstcin basis of the Miintz space f(2,2)(3) associated with the par- 
tition A — (2, 2) over an interval [a, fo], i.e., the space f(2.2)(3) = span{l, t, t^, t^). 
From Theorem [51 Lemma [1] and the branching rule (j51|) . we have 



Bli2,2){t) = y SrriU,a)SrriUM ^^^'^^^ 



where U = {a? ^ ^b^) and the Schur and the skew Schur functions can be com- 
puted, for instance, using the combinatorial definitions. 

Elementary Miintz spaces: Let (iJ^ ^^s, i?"^^^ ..., B^ ^^J be the Chebyshev- 
Bernstein basis of the rth elementary Miintz space £{it) (n) over an interval [a, b]. 
We have 

Moreover, by Lemma [TJ we have 

fx{n+l) _ /(i.)(7i+l) _n + l 
fx(o){n) fiir-i)in) r 

Therefore, a direct application of Theorem [SJ shows 

Corollary 2. The Chehyshev- Bernstein basis (_Bq ^^x, i?"(.j^rx, ..., B" .j^^s) of the 
rth elementary Miintz space over an interval [a, b] is given by 

^fe,(i'OW ^ ^k[t) e^(a"+l^^6'=)e^(a"-^6fc+l) 

Complete Miintz spaces: Let (B"^ ,^-., B^ ,..,..., B^ ,..) be the Chebyshev- 
Bernstein basis of the rth complete Miintz space £(r)(n) over an interval [a, b]. 
The branching rule (P7)) leads to 

f 5(,)(a"-^6^ ^) = ^(afo)'^-■'"i■'"/^,(a"-^fe'=). 

Therefore, applying Theorem [5] lead to the same result as in [ll[, namely 

Corollary 3. The Chebyshev- Bernstein basis (^o (r)' ^"(r)' ■•■' ^n (rl) of the 
rth complete Miintz space over an interval [a, b] is given by 

B- (f) f'' + '\l3-(t) ^UMy''h,{a-\b^)P 
'^■M^^ V n y '=^^/^,(a"+l-^6fe)/^,(a"-^6fe+l)■ 

Hook Miintz spaces: Let (-Bq mi^-,, Bw;!^-,, ..., i?^' (-n^O be the the Chebyshev- 
Bernstein basis of the {l\r) hook Miintz space £mj.\{n) over an interval [a, 6]. 
Noticing that for j = 1, .., / -I- 1, {l\r)/{j) consist of two connected components. 
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Therefore, for j — 1, ..., n, we have Sniryu-) = hiJ^-i^jCr. j = 1, ..., Ai. Therefore, 
using the branching rule (1^^ . we have 



t b(i\r) \a ,6 , — I = 6(/|r)(a ,0 )t + 



i+i 
J2hi+i-,K-\b'')er{a''-\b'^)iab)H'+'-^. 

Therefore, Theorem [5] gives 

Corollary 4. The Chebyshev-Bernstein basis of the (l\r)-hook Miintz space over 
an interval [a, b] is given by 

kAMr)\ ) ''^V + Z + IV n+1 

5(^|,)(a"+l-^6fc)5(^|,)(a"-^&fe+l) 

Staircase Miintz spaces: If iB'^^^B'^^, ..., i?" y) is the Chebyshev-Bernstein 
basis over an interval [a, 6] of the ^-staircase Miintz space E\(ri)^ where the 
partition A is given in ((39]) . then from ((38|) . we have 

5?,AW = i3^'(t'+'), for fc = 0,...,n, 

where i3[? is the classical Bernstein basis over the interval [a'+^,&'+^]. Our 
objective here is then to show that Theorem [5] reconfirm this fact. The method 
of computation consists in using equations (|47p for the partitions A and A*-*^' 
and inserting these equations into Theorem [SJ We will omit all the details of 
the computation but only mention that it is helpful to rename (a"^'^,6'"') — 
(ui, ...,Un) in order to detect easily the simplifications and that by equations 
(|T7)) . the term 

/,<o,(n) 5„o,(l") ^ ^ ^ 
appears naturally within the computations. We find 

_ fn\ it - fl)^-(6 - t)»~fe hi{t,a)^^hi{t, br-'^ 
''■'^^ ' \kj (fe-a)" hi{a,b)^ 

Therefore, 

„ /n\ (t'+i-a'+i)fc(&'+i^t'+i)"-fc 

as expected. 

7.1. Weighted de Casteljau paths 

Our strategy for computing the Chebyshev-Bernstein basis associated with 
a Miintz space £\{n) consists of computing the product of weights associated 
with specific paths in the de Casteljau algorithm. Working directly with the 
Pascal-like graph of the de Casteljau algorithm reveal to be challenging and 
induction arguments does not seem to work. For these reasons, we will define 
a combinatorial object that, in some sense, can be viewed as one-dimensional 
projection of the two-dimensional paths in the de Casteljau graph. 
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Definition 2. A set A„ = (Aq, Ai, A2, ...., A„) is said to be a de Casteljau path 
in {0, 1, ..., n} if and only if, 

i) Every set Ai is a subset of {0, 1, ..., n} such that \Ai\ = ^ + 1. 

ii) The set Ai^i is obtained from Ai by adding to Ai an element of the form 
max{Ai) + 1 or min{Ai) — 1 under the condition that the set remains a subset 
o/{0,l,-,n}. 

We will often represent a de Casteljau path as 

^0 >Ai > An-l > An 

in which the subsets Ai,i = 0, ...,n are viewed as vertices and the arrows are 
viewed as edges. We will also adopt the convention of writing the elements of 
the set Ai,i = 0, ...,n in increasing order. Note that for any de Casteljau path 
A„ = (Aq, ..., An) in {0, 1, ..., n}, we necessarily have An — {0, 1, ..., n}. 

Example 4. Examples of de Casteljau paths in {0,1,2,3} and {0,1,2,3,4} 
respectively can be given as 

{1} -^ {0, 1} -^ {0, 1, 2} -^ {0, 1, 2, 3} 



or 



{2} -^ {1,2} -^ {1,2,3} -^ {1,2,3,4} -^ {0,1,2,3,4} 



Definition 3. Let a,b and t be real parameters and let ■0 = (V'^,V' ) be a 
non-zero two components real function 

V'(ui,...,u„-i;a, 6,i) = (V'+(wi,...,w„-i;a, &,i),V'~(wi,...,w„^i;a, &,i)) 

where ■0+ and ip~ are symmetric in the variables Ui,i — l,...,n — 1. A ip- 
weighted de Casteljau path A„ = {Aq, Ai, ..., An) is defined as associating a 
weight to every edge Ai — > Ai^i of the path A„ according to the following rules 

1) If Ai-^-i is obtained from Ai by adding the element max{Ai) + 1, then the 
weight of the edge is given by 

V'+(tl^'l^\ fe^^f-^'+i^ a"-l^'l-™*"(^'+i); a, b, t). (52) 

2) If Ai-f-i is obtained from Ai by adding the element min(Ai) — 1, then the 
weight of the edge is given by 

^-(^|A,|-l^ ^™n(A, + i)^ a"-|^<|-™m(A, + l). a, b, t). (53) 

We will represent the weight on the edges as 

Ai^Ai+, or Ai^Ai+, 

in which there is no need to write the arguments of the functions ip^ and ip~ as 
they are uniquely defined from the sets Ai and Ai^i according to the rules \5'^] 
and {23). 

We define the weight W^_a„ (a, b, t) of a ip-weighted de Casteljau path A„ as 
the product of the weights of the edges. 
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Note that by a simple application of the pigeonhole principle, for every de 
Casteljau path A„ — {Aq, Ai, ..., Ai), we have \Ai\ + min{Ai+i) < n. Therefore, 
all the exponents (referring to the number of occurrence of the arguments) in 
((5^ and ([55)1 are non- negatives. 



us define a function ■01 = [ij)^ , ■0]^ ) by 



Let A be a partition of length at most n and A'*'' its bottom partition. Let 

ion ■01 = {ip^,%})^) by 

ip'^{ui,...,Un-.i\a,b,t) 



l + (.. .. _ i_ _^^ _Sx(a){ui,...,Un-l,a) 



S\(0){ui,...,Un^l,t) 

V'l (■ui,...,u„_i;a,&,t) =--^ -r. 



We have 



Proposition 9. The weight W^^^^ (a, fe, i) o/ any de Casteljau path A = (^o, •■•, An) 
with Aq — {k} and ipi defined in |5^[ ) is given by 

W (n h i\ S^(o){a''^'',b'') 
l4^A,^i (a, 0, t) = J — T . 

Proof. Let A„ = {Aq,Ai, ..., An) be a de Casteljau path and consider the generic 
product of weights over two arbitrary adjacent edges 

Ai — > Ai+i — >Ai+i. 

In general, we would have four situations with regards to the weights, namely 

At ^ At+, ^ At+,, Ai ^ Ai+, ^ Ai+2, (55) 

Ai ^ A,+, ^ Ai+2, Ai^A,+,^Ai+„ (56) 

Let us consider the first case Ai —^ ^;+i —^ A1+2 and let D be the denomi- 
nator of the first edge and N the numerator of the second edge. By definition 
131 we have 

D = S'A(o)(il'^'l~\6"""(^'+i\a"~l'^'l""""('^'+i\i) 

and 

In this case, we have min{Ai) — min{Ai-^-i) = min{Ai^2)- Counting the number 
of occurrence of the variables t, a and b in the expression of D and N, shows that 
D — N . Similar arguments show that also for the other situations in (j55p and 
([5^ the denominator of the first edge is equal to the numerator of the second 
edge. Therefore, upon taking the product of the weights of a de Casteljau path, 
the denominator of the weight of an edge will be simplified with the numerator of 
the weight of the adjacent edge. The two factors that survive the simplifications 
are: the numerator of the weight of the first edge and the denominator of the 
weight of the last edge of the de Casteljau path. Let us compute the numerator 
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of the weight of the first edge: Since, by assumption, we have Aq — {fc}, we will 
have two situations 

Ao = {k}^{k,k + l} or Ao = {k}^{k,k-l}. 

In the first case, the numerator is Sx(a){t'~',b'',a'"'^'^^^,a), while in the second 
case, the numerator is 5;^(o) (i°, 6^^"^, a"~'^,6). Therefore, in both cases, the 
numerator is Sx{o){a^~'',b'^). For the denominator of the weight of the last 
edge, we can only have a single situation, which is 

A„-i^{0,l,...n}, 

and in which the denominator is given by Sxm (i")- D 

Remark 6. Note that we did not use the fact that S^m is a Schur function 
and all what was needed is for the function S';^(o) to be a symmetric function. A 
similar remark can be applied to the next proposition. The reason for us not to 
state the propositions in their full generality is to make the exposition for later 
use more transparent. 

Let A be a partition of length at most n and consider the function il)2 = 
{tp 2 ,1^2) defined by 



V'^(ui,...,u„_i;a, 6,t) 



Sx(ui,...,u„-i,b,t) 



Sx{ui,...,Un-i,a,b) 

V-z (ui,...,u„_i;a,6,t) =-— -. 

i)A(wi,---,Wn-i,a, 0) 

In this case, there would be simplifications in the weight associated to every de 
Casteljau path, however there is no simple close explicit formulas that embodies 
all the simplifications. In the following, we will show that the specialization t — a 
or i = 6 on the weights of every de Casteljau path is given by a simple closed 
formulas. More precisely, we have 

Proposition 10. For any de Casteljau path A„ = (Aq, Ai, ..., An) such that 
Aq — {k}, we have 

Wi^^A.)i^^bM=a= g^^^l.J^,^ (58) 

and 

W,^,,^^,iaAt)U=,^ ^^^^,l,J^,y (59) 

where -02 is defined in ( [57p . 

Proof. Let us start with the first equation (j58p of the Proposition. Noticing 
that for t = a, we have ip2 = 1 shows that the edges with weight 1 does not 
contribute to the total weight of a de Casteljau path. Let A„ — {Ao,Ai, ..., An) 
be a de Casteljau path and consider a situation in which a part of the path has 
the following weights 

Ai ^ Ai+i 4 Ai+2 4 .... 4 Ai+h-i ^ Ai+h (60) 
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with h > 1. Consider the numerator N of the weight of the first edge of (pO)) 
and the denominator D of the weight of last edge of (|M)) . We have 

AT = 5-^^(^1^11-1, &"""(^'+i),a"-l'4il~™"(^i+i)^a^i) 
and 

From the structure of the path (|60p . we can see that min{Ai^fi) = min{Ai^i) — l. 
Therefore, the number of occurence of 6 in iV is equal to the number of occurence 
of fo in _D. This shows that when t = a, we have N — D. Thus, for any de 
Casteljau path, the numerator of an edge with weight 'ip2 '^ill be simplified 
with the denominator of the weight of the closest edge with weight "4^2 ■ There 
is two special cases that will need a separate treatment. The case in which there 
is no edge with weight 1^2 ^^'^ the case where there is a single edge with weight 
ip2 ■ In the former case, we only have a single de Casteljau path, namely 

{0} ^ {0, 1} 4 {0, 1, 2} ^ .... 4 {0, 1, 2, ..., n} 

with weight equal to 1. This case corresponds to a de Casteljau path A„ — 
(^0,^1, ■■■,An) with Aq = {0} and the weight of the path is consistent with 
equation (1551) of the proposition for fc = 0. In the latter case, there exists an 
r < n such that the de Casteljau path will look like 

{1} ^ {l,2}....{l,2,...,r} ^ {0,l,...,r} ^ .... 4 {0, 1, 2, ..., n}. 
In this case the weight of the path, at the specialization i = a, is given by 

As this case corresponds to a de Casteljau path A„ = (Ag, Ai, ..., A„) with 
Aq = {1}, the result is again consistent with the first equation of the proposition 
when k — 1. For all the other cases, the weight of the de Casteljau path is then 
a ratio in which the denominator is given by the denominator of the first edge 
with weight '0^ , and the numerator is the numerator of the last edge with weight 
ip2 ■ For the first edge we can only have the following two situations 

{fc}^{fc,fc-l} 
and 

{k} ^ {k, k+l} ^ {k, fc+1, fc+2}.... \ {k, fc+1, fc+2, ...} ^ {fc-1, fc, fc+1, ...}. 
In both cases, the denominator, when i = a, is given by 

For the last edge, we would have the following situation 

Ai ^ Ai+i ^ Ai+2 ^ .... 4 {0, 1, ..., n}. 
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In this case, we have min{Ai^i) = and the numerator of the edge Ai — ^ ^i+i, 
when t = a, is given by 

^A(a"+'). 

This prove the statement of (|58p. Similar arguments when we take the special- 
ization t = b lead to the second equation ([5^ of the proposition. D 

Consider, now, the triangular Pascal-like graph of the de Casteljau algo- 
rithm. We encode every node of the graph as follows : the node that lies in 
the rth horizontal level and the fcth position going from left to right is encoded 
as {k — l,k,k + 1, ..., k + r — 2}. For example, the code associated with the de 
Casteljau algorithm based on 4 control points is given by 

{0} {1} {2} {3} 

\ /\ /\ / 

{0,1} {1,2} {2,3} (r) 



{0,1,2} {1,2,3} 

\ / 

{0,1,2,3} 

Let A be a partition of length at most n and denote by ■0 = [tp^ ,ip~) the 
function such that 

?/;+(ui,....,w„_i;a, &,i) =^(ui,....,u„_i;a, 6,i) 

(61) 
ip (ui,....,-u„_i;a, 6,i) =a(ui, ...., u„_i; a, 6,i), 

where a and /? are the pseudo-afhnity factors of the space £\{n) as defined in 
Theorem 31 We have 

^ = (V'+, v") = {^i^t^t, ^V'rV'2"), (62) 



where f/'i and '(/'2 are defined in ((54)) and ([57| respectively. Consider the de 
Casteljau algorithm based on the control points ((^9(0"), 1^9(0""^, 6), ...., (y9(fo")) 
where ip is the Chebyshev blossom of a Chebyshev function (f> over an inter- 
val [a, &]. Let us write a generic triangle in the de Casteljau algorithm with 
the value of the pseudo- affinity on the edges of the triangle and write the 
same triangle with our coding of the de Casteljau algorithm and in which the 
weight on the edges follow the rules ([5^ and ([55]) of Definition [3] We have 

(^(a"-'-'=+2,&'='i,r-i) ¥3(a"-'-'=+\ 6*^,^-1) 
pf^^n-r^k+i^ fe'=-i , r-i\a, h, t) a(a"-'-y 1 ,h^-^ , T^i ; a, 6, t) 

-r-fc+l hk-1 



and 
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{k- l,k,...,k + r -2} {k,k + 1, ...,k + r - 1} 
V'+(a"-'-'=+\ fe*^-!, rV; a, b, t) V'" (a'^-'^V+i, fo*^"!, T^i; a, 6, t) 

{k- l,k,...,k + r - 1} 



Realizing that the weights in both of the edges of the triangles are the same 
shows that if we denote by (Sg A' ^i A' •■•' ^n a) ^'^^ Chebyshev-Bernstein basis 
of the Miintz space £\{n) over the interval [a, 6], then the de Casteljau algorithm 
is the claim that for A; = 0, ..., n 

A„ 

where the sum is over all the de Casteljau paths A„ — {Aq, ...,An) such that 
Aq = {k}. It is simple to see from ((62)) that for any de Casteljau path A„ — 
{Aq, ..., A„) such that Ao — {k}, we have 

(b - tY^'^it ~ aY 

W^V,A„(a,^,^) = ., _ ^„ WVi,A„(a,6,t)Wv,2,A„(a,&,i) 

[0 a) 

Using Proposition ini for W^j.A„(a,&, i), we obtain 

Proposition 11. Let X be a partition of length at most n and let (i?o ;^, B^ ^^, ..., 
^n \) b^ ^^^ Chebyshev-Bernstein basis of the Miintz space £\{n) over an inter- 
val [a, b]. Then, we have 

k,\\ J [b-aY S)^(o){t") ^^ *'' "^ ' ' '' 

where the sum is over all the de Castelaju paths A„ = (Ao,...,A„) such that 
Aq — {k}. A*^"' is the bottom partition of the partition A and the function ip2 is 
defined in ( |57| ). 

There is two special cases in which the quantity W^^^Aio,,b,t) can be given 
an explicit expression. Let us consider the set of the de Casteljau paths A„ — 
{Aq, ..., An) such that Aq ~ {0}. In fact, there is a single path which is 

{0} ^ {0, 1} ^ {0, 1, 2} ^ .... ^ {0, 1, 2, ..., n}. 

As the reader can readily check, the product of weights along this single path 
is given by 

For similar reasons, there is a single path A„ = (Aq, ..., An) such that Aq = {n}, 
namely 

|n| > |n — l,n| > {n ~ 2,n — l,n} 5- .... > |L), 1, 2, ..., n|. 
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ib- 


-t)" Sx{b,n S^m{a^) 


(b- 


-a)"5A(&,a")^A(«)(<") 


(t- 


- a)" 5a (a, t") 5^(0,(6") 



Along this path we have 

. , N -5'A(a,r) 
5a (a, 6") 

Therefore, according to Proposition II 11 we have 

Proposition 12. Let X be a partition of length at most n and let (Sg aj ^"a' •■•> 
S^j j!^) be the Chehyshev-Bernstein basis of the Miintz space £x{n) over an inter- 
val [a,b]. Then, we have 

Blxit) 

and 

"^^^^ (6-a)»5A(a,&«)5A(o,(i")- 

In general, the expression of the Chebyshev-Bernstein basis obtained by 
computing the weight on the de Casteljau paths has complicated expressions. 
Let us for example consider a partition A of length at most 2 with it associ- 
ated Chebyshev space £\{'2). Proposition [T^] provides us with the Chebyshev- 
Bernstein functions B^ ^ and -B| a • In order to compute Bf ^, we should com- 
pute ^j^ W^2_A2(aj b, t) where the sum is over all the de Casteljau paths A„ = 
(^0,^1,^2) with ^0 = {!}■ In this case, we have two de Casteljau paths 
namely, 

{1} ^ {0,1} ^ {0,1,2}, {1} ^ {1,2} ^ {0,1,2}. 

Computing the weights along these two paths lead to 
2 u\ ib-t){t~ a) Sx(o){a,b) f Sx{a,a,t)Sx{t,t,b) Sx{b,b,t)Sx{t,t,a) 



[b-aY Sxm{t,t) \Sx{a,b,t)Sx{a,a,b) Sxia,b,t)Sx{b,b,a) 

It is rather surprising that with this expression in hand the function Bf a is a 
polynomial function in t. 

Although summing the weights over the de Castelaju paths does not lead to a 
practical method of computing the Chebyshev-Bernstein basis, the concept leads 
to the following important information about the derivatives of the Chebyshev- 
Bernstein basis, which by some inductive argument on nested Miintz spaces will 
provide us with the desired explicit expression. 

Theorem 6. Let A = (Ai, A2, ..., A„) be a partition of length at most n and 
let {BJ} x{t),B^ xi^)' ■■■T^n \i^)) be the Chebyshev-Bernstein basis of the Miintz 
space £\{n) over an interval [a,b]. Then, we have 

(,) n\ a^^ /A(n + l)5A(c)(a"-^&'-) 
Bk,x (a) 



and 



{n-k)l{b-aY /a(o)H 5A(a"+l-^fe^) 
S?A^""'^^W = (-1) 



("-'^)m - . .^n-fc "! b^^ /A(n-H)5A(o)(a"-^6'') 
kl{b-ay^-^^ fxmin) 5a (a"-^ &'=+!) ' 

where X^*^' is the bottom partition of X. 
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Proof. Let us define the function H{a^ 6, t) by 

C f^n — k i^k\ 

H{a,b,t) = ''^'°'^° ,^X ^ E^^^.A„(«,&,t), 

where the sum is over all the de Casteljau paths A„ — {Ao,Ai, ..., An) such that 
Aq — {k}. From Proposition [TTJ we have 

^k,x{i) JJ^ZTay^ H[a,b,t). 

By the Leibniz derivative formula, we then have 

Bl>}''\a) = -j^ "^ -HiaAth=a. 

(J(n-fc)!(6-a)'^ 

Using Proposition [TUl the fact that there is C^ de Casteljau paths A„ — 
{Aq, ..., An) such that Aq = {k} and the fact that for any partition jj,, we have 
Sfj,{x"') = ffj,{n)x^'^^, we arrive at the first equation of the Proposition. Similar 
treatment on the {n — k) derivative of the Chebyshev-Bernstein basis at the 
parameter b lead to the second equation of the proposition D 

7.2. A Descent Construction of the Chebyshev-Bernstein Basis 

In the following, we will show that Theorem [5] will allow us to relate the 
Chebyshev-Bernstein bases of Miintz spaces associated with two different parti- 
tions under a condition on the partitions that we state in the following definition 

Definition 4. Let X be a partition of length at most n. A partition fi of length 
at most (n + 1) is said to be a dimension elevation partition of X if, and only if 
the Miintz spaces associated with the partitions X and fi satisfy 

£x{n) (Z£^{n + l). 

We can characterize all the dimension elevation partitions of a given partition 
A as follows : 

Proposition 13. Let X = (Ai, ..., A„) be a partition of length at most n. Then 
every dimension elevation partition pi is of the form 

/i = (r + Ai,r + A2, ...,r + A„,r), with r > 0, (63) 

or of the form 

M = (Ai - I, A2 - 1, ..., As - l,p, As+i, ..., A„) with 1 < s < n, (64) 

under the condition that pL is a partition. 

Proof. Let (j) be the Chebyshev function associated with the partition A. Then 
4) is given by 



m = (i^^ 



-A2 + I .A1-A3+2 ,Ai-A„ + (n-l) ,Ai + 
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There is two ways to supplement the function cf) with an extra component of the 
form t™ with m a positive integer. 

The first way: We can add a function i™ such that m is strictly larger than 
all the exponents in the components of the function (f), i.e., ttt, > Ai + n. In this 
case, if we denote by /i = (^i,/i2, ...,/^n+i) the partition associated with the 
obtained space, then we have 

m = ^ll + {n+l), Xi-Xi = ^ii-^ii, for z = 2, ..., n, and ^i-/i„+i = Ai. (65) 

If we denote by r = /ii — Ai, then r > as tti > Ai + n. Moreover, equation 
(|65p shows that we can write ^^ as r + A^ for i = 1, ..., n and fin+i — r, thereby 
leading to the form of the partition in (j63p . 

The second way: We can add a function of the form t™ in which the exponent 
m lies between two exponents of the components of the function </), i.e., for an 
1 < s < n, we insert t™ between t^i-^s+(s-i) g^^^^ ^Ai-a,+i+s (^Xn+i = 0). 
Note that this is possible only if A^ > A^+i. In this case, we would have 
Ai + n = /ii + (n + 1) and therefore, /ii = Ai — 1. By using the condition that 
we should have 

Ai-Aj=^i-^i, for z = 2,...,s 

and 

Ai - Ai - 1 = ^1 -^i+i, for i = s + l,...,n. 

we arrive at a partition of the form ([M]) with p = Ai + s — (to + 1). D 

Consider, now, a partition A of length at most n, and let /i a dimension 
elevation partition of A. As an element of £fi{n + 1), the Chebyshev-Bernstein 
basis BJ^ ^, k — 0, ...,n, oi £\{n) over an interval [a, 6] can be expressed as linear 
combination of the the Chebyshev-Bernstein basis B^^^,k — 0,...,n + l, of 
£fj.(ji + 1) over the same interval. Exhibiting the vanishing properties of the 
Chebyshev-Bernstein bases as expressed in Theorem [T] shows that [2| 

As Proposition [5] gives explicit expressions of all the needed derivatives in the 
last equation, we can express the Chebyshev-Bernstein basis associated with 
the partition A in term of the Chebyshev-Bernstein basis associated with a 
dimension elevation partition [i,. To write the expression in a more compact 
fashion we define the following factors 

Definition 5. Let A (resp. fi) be a partition of length at most n (resp. at most 
(n + l)). Denote by X^^' (resp. ^^^') the bottom partition of X (resp. ^). For 
k = 0, ..., n, we define the following factors 

and 

A-in k) ^^'°' ("""'' ^')^^ (""''' ^'^') (68) 
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(69) 



A A hm{n)Uin + 2) 



Now using equation (j66p and in which we insert the value of the needed 
derivatives from Proposition [6] we arrive at the following 

Theorem 7. Let X ~ (Ai, ..., A„) he a partition of length at most n and let /x = 
{^Li, ..., ^n+i) be o, dimension elevation partition of X. Denote by {B^ ^, ..., B^ ^) 
(resp (-Bqu !---i^n+iu)J ^^^ Chebyshev- Bernstein basis of £\{n)(resp. f^(n + 
\)) over an interval [a, 6]. Then, we have 

where p — Xi — fii. 

To illustrate the use of the last Theorem as a mean of finding explicit ex- 
pression for the Chebyshev-Bernstein basis, consider the rth elementary Miintz 
space i.e., the Miintz space associated with the partition A = (f^). The parti- 
tion /i — (0) is a dimension elevation partition of A whose Chebyshev-Bernstein 
basis over an interval [a, b] is given by the classical polynomial Bernstein basis 
of order n + 1 over the interval [a,b]. We have A*^°^ = {V~^), /i = /xf^"^ = (0) 
and ( ) = (n + l)/r. Therefore, applying Theorem [7] leads to 

_ e.-l(a"-^6'') / (n + 1 - k)a +i {k + 1)6 „„+i,,A 

Bkxm W - - [-^^^TT^k;b^)^k ^*) + ^;(^^^^k;b^^^+^ ^^>) ■ 

This illustrative example prompt us to consider the following algorithm for 
computing the Chebyshev-Bernstein basis associated with a partition A. We 
can construct a sequence of nested Miintz spaces ^p^^., (n -I- j),j = 0, ..., m such 
that 

£\{n) = ^p(n)(n) C f^(j)(n + l) C .... C fp(^^_j,(n+TO-l) C £:^(^,(n+TO), (70) 

where the partition /i(m) — 0. As the Chebyshev-Bernstein basis of the space 
£fj,,s (n + m) is the classical Bernstein basis of degree n + m, we can construct 
iteratively the Chebyshev-Bernstein bases starting from the space £^ ^^_ (n -|- 
m — 1) until reaching the space £\(n) using Theorem [71 

There is several sequences of nested spaces that satisfies (jTO)) . starting from 
the space £\{n) and in accordance with the constraints of Proposition [T21 For 
example we have the following two sequences of nested Miintz spaces 

^rm(»^) C £^n + 1) C £'n(n -I- 2) C £0(n -f 3) 

and 

g rrn (?^) C £rHn + 1) C fn(" + 2) C f0(n + 3) 

which correspondent respectively to the following situation of nested Miintz 
spaces 

(1, i^, i^ t\ ..., i"+3) c span{i, i^, i^ i^^ i^ ..., ^"+3) c 

span(l, t^ t^, f*, t^, t^ ..., r+3) c span{\, t, t^, t^ t^ i^ i^ ..., r+3) 
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and 

(1, t^, t\ f, ..., r+3) C span{l,t, t^, t^,t^, ..., e+^) C 

span{i,t, t^ i^ i^ i^ ..., r+3) c span{i, t, t^ , t^ t^ t^ t^ ..., r+3). 

Remark 7. In some circumstances, it is not necessary to have a full descent 
of nested Miintz spaces as in ([70)) to compute the Chebyshev-Bernstein basis 
of a specific Miintz space. We can sometimes use the staircase Miintz spaces 
as a short-cut space for the computation. For example, consider the Miintz 
space E ~ span{l,t^,t^,t^). The partition associated with this space is given 
by A = (5, 3, 1). To compute its associated Bernstein-Chebyshev basis, we can 
make the dimension elevation E = span{l,t'^,t^,t^) C F = span{l,t'^,t'^,t^,t^). 
As the Chebyshev-Bernstein basis of F over an interval [a, b] is known in terms 
the classical Bernstein basis over the interval [a^,6^], we can use Theorem [7| 
to find the Chebyshev-Bernstein basis of the space F in a single iteration. In 
principle, we can use this trick to compute the Chebyshev-Bernstein basis of 
any Miintz space whose components are a "reparametrization" , in the sense of 
((55)l . of an already studied Miintz space. 

In the following, we will show that there is a particularly convenient choice 
of a sequence of nested Miintz spaces in which an inductive argument along the 
sequence will give us the explicit expression of the Chebyshev-Bernstein basis. 

Definition 6. Let X be a partition of length at most n. We define the border 
complement rj of the partition A as the partition obtained by removing the first 
column of X. In other word, if X is given by X = (Ai, A2, ..., A^, 0, ..0) where 
Ai > 1 for i — 1, ...,s, then its border complement is given by rj = (Ai — 1, A2 — 
1,...,A3-1,0,...,0). 

It is clear from proposition [T3] that if A is a partition of length at most n 
and 77 its border complement then 77 is a dimension elevation partition of A i.e., 
£x{n) C £r}{n + 1). A first hint of the usefulness of this choice is the following 
combinatorial lemma 

Lemma 2. Let X be a non-empty partition of length at most n and let rj be it 
border complement. Then, we have 

V\ .f\{n + l)f^(o){n + l) n+1 



vA;„ ^(n + 2)/;,(a,(n) /ia(1,1)' 

where h\{l,l) is the content of the first square of the partition X. 
Proof. From Lemma [U we have 

hjn + l) _-^ (n + l) + CA(l,j) 
/a(o)H f\ hxil,j) 

and 

fvin + 2) ^ ^^' (n + 2) + c,(l,j) 
^(o)(n + l) ijL h^ilj) 

Moreover, from the definition of 77, we have for j = 1, .., Ai — 1 

c„(1,j) = ca(1,J + 1)-1 and h^{l,j)^hx{l,j + i 
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Using this extra information in the computation leads to a proof of the lemma. 

D 

With a border complement partition as a dimension elevation partition, The- 
orem [7] takes the simpler form 

Proposition 14. Let X be a partition of length at most n, and let rj be its border 
complement. Then, the Chebyshev- Bernstein basis associated with X and rj over 
an interval [a, b] are related by 

We need the following proposition and the subsequent corollary in order to 
give a proof of the main Theorem [5] 

Proposition 15. Let X be a partition of length at most n and fx it border 
complement. Then, for any real numbers U — (ui, ...,u„„i), and real numbers 
X and y, we have 

Sf,(o){U,x)Sx{U,y) - S^(o){U,y)Sx{U,x) = {y - x)S^,{U,x,y)S)^(o){U), 

where X^^' (resp. /i'^V the bottom partition of X (resp. fi). 

Proof. Let us first assume that the partition A is of exact length n. Using ([T7)). 
we have 

Sx{U,y) = y(Y[uA S^{U,y), SxiU,x) ^x(Y[uA S,,{U,x), 

and 

5a(o)((7)- \J{uAs^,.,{U). 

In this case the statement of the proposition is nothing by Proposition [71 Now, 
let us assume that A is of exact length k < n i.e., A = (Ai, ..., Afc,0, ..., 0) with 
Afc > 1. Consider the Miintz tableau (/x'^''\ /x'^-'^^ ..., /x^"^) associated with the 
partition /i. Then, we have, fi^'^' = A. Therefore, applying Proposition [5] to the 
partition fi leads to a proof of the proposition. D 

Corollary 5. Let X — (Ai, A2, ..., A„) be a partition of length at most n and fi 
it border complement. Then, for any positive real numbers U = (mi, ...,u„_i), 
and positive real numbers a, b and t, we have 

a(b - t)t^'-^S^{U, a, ^)Sx{U, b) + b{t - a)t^'-^S^{U, b, j)Sx{U, a) = 

{b~a)t'''S^{U,a,b)Sx{U,^). 

Proof. Applying Proposition [T51 with x = a and y = ab/t and multiplying both 
sides by t^^ give 

a{b - t)t^'-'^S^,{U, a, ab/t)S^m (U) = 

t^' {S^^oAU,a)Sx{U,ab/t) ~ S^,o>{U,ab/t)Sx{U,a)) . 
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Therefore, we have 

a{b -t)t^'-^Sf,{U,a,ab/t)Sx{U,b) = 

^Ai MIA (^ ([/, a)Sx{U, ab/t) - S^m {U, ab/t)S,{U, a)) . (71) 

Similarly, by applying Proposition 1151 with x — b and y — ab/t, we have 
b{t ^ a)t^^-'^S^{U,b,ab/t)Sx{U,a) = 

*^' M^ (5 <o, ([/, a6/i)5A(t/, fe) - 5^(0, (t/, fe)5A([/, aVO) • (72) 

Summing up the expressions ([7T]) and ([7^ shows that the right hand side of the 
equation in the corollary is given by 

^^'?^^Vn^^^ (^^('^' ^)^M<») iU, a) - Sx{U, a)5^,o, {U, b)) . 
Oxm[U ) 

Using again Proposition [151 to the quantity between the parenthesis in the last 
equation, with x = a and y — b, leads to the desired formulas D 

Proof of the main Theorem [5} 

Proof. We will proceed by induction on the number of boxes in the partition A. 
For an empty partition, the Chebyshev-Bernstein basis is given by the classical 
Bernstein basis which is consistent with the formula given by the Theorem. 
Let us assume that the Theorem is true for any partition with less than m 
boxes. For a given partition A with m boxes, let us denote by /i its border 
complement partition. Then, necessary the number of boxes in /i is less than 
m. By Proposition 1141 we have 

{n + l-k)a 5,(o,(a"-^&'=)^>;^+2-^6^ 



^k,\\^) - T:-TT~r\ C ^ l„n+l-k hk\C.I„n+l-k l,k\^k.M ^^1^ 

(73) 



hx(l,l) S^^o){a^+^-\b^)Sx{d-+^-\b^) "■■'' 



(fc+i)6 5A(o)(a"'^fe^)^p(«"-^&''+^) p„+i ,.. 
hx{i,i) v(«""^''''+')'^A(a"-^&'=+l) fc+i,M^ )■ 

Using the induction hypothesis (l48l) on the Chebyshev-Bernstein functions B'^^l {t) 
and B^^^{t) and carrying out all the obvious simplifications as well as using 
Lemma [21 we find that the first term in (j73p is given by 

fxjn + 1) a{b ~ t) „ . . A,-i ^A(o) (a"-^ &")5^(a"+l-^ b\ ab/t) 
fxmin) b-a ''^ ' 5A(a»+l-^ 6'=)5^(a"+l-^ fe'^+i) ' ^ ' 

while the second term in (j73p is given by 

fxjn + 1) bjt - a) „. . A,-i ^A(o) (a"-^ 6")5^(a"-^ b^+\ab/t) 
/A(o,(n) b-a ^^^ ^A(a"-^6'=+l)5^(a"+l-^6''-+l) ' ^ ' 

Summing the two last equations leads to 

~ /x\n) " ^^^M(a"+l-^6^•+l)^^ ' ' >' ^^^' 
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where il is given by 

f t^'-\a{b-t)S^{a^+^-\b\ah/t)Sx{a^~\h^+^) 

lX{a,0,t) ^_^[ ^;,(a»+l-^6fc)^^(a»-^6fe+l) 

b{t - a)5'^(a"-^ &^+\ aVt)5A(a"+^~^ 6'') 

Using Corollary [SI with U ~ (a"^''', 6*^) (note that here we view A as a partition 
of length at most (n + 1) to be able to take U with n components) shows that 

Inserting the last term into equation (l76l) result in the proof of the Theorem. D 

Remark 8. One we have guessed the explicit expression of the Chebyshev- 
Bernstein basis, it is in principle, possible to find a simpler proof than the one 
given here. For instance, according to the characterization of the Chebyshev- 
Bernstein basis given in [9| and in view of Proposition [SI we need only to show 
that every element B^ -^ expressed in (PSJ) is an element of the Miintz space 
£x(n). However, the advantage of our proof lies in demonstrating the elegant 
combinatorics beneath the relations of Chebyshev-Bernstein bases associated 
with different partitions. 

7.3. Dimension elevation process 

Let A — (Ai, A2, ..., A„) be a partition of length at most n and let 77 = 
(771,772, ■•■, '7n+i) be a dimension elevation partition of A. Consider a £\{n)- 
function P written in the Chebyshev-Bernstein bases associated with the par- 
titions A and 77 over an interval [a, b] as 

n n+1 

P{t) = ^ Blx(t)Pu = Y. BZ'(t)P^- (78) 

fc=0 k=Q 

Using Theorem [7] to detect the coefficients of B^'\^{t) in the expansion (l78l) . we 
readily find 

Theorem 8. The Chebyshev-Bezier points Pk in |y<$[ ) are related to the Chebyshev- 
Bezier points Pk by the relations 

Pq = Pq, Pn+l = P71, 

and for k — 1,2, ...,n 

Pk == P[\,r)]{n, fc - 1) Pk-i + ^[x,rj]{n, k) Pk, (79) 

where C[a, ?;](", k) and p^\.rj\(n, k) are given by 

(ri + l-fc)a^^ /r7\ 
C[x^^]{n,k) = —^ l^^j ^l{n,k) 

and 

P[x,r,]{n,k) = -^^ i^^j Al{n,k), 

where r^(n, fc) and A^{n,k) are defined in |6'7| j and i68\) . 
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Figure 1: The dimension elevation process £ri\{3) = span(l, t'^, t^, t*) C ^0(4) = 
span(l, t, t^, t^, t^). (Poi ^i> ^2, ^3) are the Chebyshev-Bezier points of a £n)(3)-function, 
while {Pot Pi, P2y P:it P4,) are the Chebyshev-Bezier points of the same function viewed as a 
^0 (4)-function. (see example 11). 



Remark 9. As the relation ([7^ is independent of the £> (n)-function P and in 
view of Q, we have px^^iin, fc — 1) + S,\^n{n, k) — 1. This relation can also be 
directly proven (with rather great efforts) using Proposition 1151 

Example 5. Let P he a £rir-^{n)-function. The partition rj ~ (0) is a dimension 
elevation partition of X. Therefore, the function P can be expressed as 

n n+1 



fc=0 



k=Q 



In this case, we have 



P[1M,(0)](", fc - 1) = -—-^—^^^-^^ 



and 



C[H-Kmin,k) 






Therefore, from Theorem [H we have Pq — Pq and Pn+i — Pn o,nd for k 
l,...,n, 

~ _ fc6e,_i(a"-'--+i,&fe-i) (r^ + l-fc)ae,_l(a"-^6fe) 

The case r = \ provide us with the following simple relationships 



Pk = 



kb p , {n + 1- k)a 



{n + l-k)a + kb (n + I - k)a + kb 

Figure[l\ shows an example of dimension elevation process for the case r = 1. 
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Pa = Pa 



Pi = P, 



Figure 2: The dimension elevation process f(i,i)(3) = span{l,t^,t'^,t^) C £'(i)(4) = 
span(l, t^, t"^, f*, t^). (Po, Pi, P2, ^3) are the Chebyshev-Bezier points of a £n j)(3)-function, 
while (Po, ^ii ^2, ^3, ^4) are the Chebyshev-Bezier points of the same function viewed as a 
fn ) (4)-function. (see example 12). 



Example 6. Let us consider the dimension elevation process 

£(i){n) C f(i_i)(n-|-l). 
Then, if we write 

n n+1 



p(t) = E Biii)it)Pk = J2 B'^AiUim 



fe=0 



fc=0 



would ha 



P[{i),{i- 



.i)](n,fc-ij-- _-___^^ 



and 



^m,(i~i)]{n,k) 



(n + l-A:)a/ii_i(a"+2-'=,&'=) 



(80) 



I hi{a-^-^+^ ,b^) ' 

Therefore, we have Pq ~ Pq and Pn+i = Pn o-nd for fc = 1, ..., n, 

Figure\^ shows an example of the dimension elevation process for the case 1 = 2. 



Pk 



8. Toward shaping with Young diagram 

A polygon V = {Po,Pi, ...,Pn) can be viewed as the control polygon of 
a £A('^)-function, where A is a partition of length at most n. Therefore, by 
varying the partition A, the curve associated with the control polygon will also 
vary accordingly. In such circumstances, the Young diagram can be viewed 
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Figure 3: The effect of iteratively adding boxes to the first row of the partition A = (2, 1). 
The black curve refers to the Chebyshev-Bezier curve over the interval [a, b] = [1, 4] associated 
with the partition A, the red curve correspond to the partition (3, 1), the green curve to the 
partition (4, 1) and the blue curve to the partition (5, 1). 





Figure 4: The effect of iteratively adding a column to the partition A = (2, 1). The black 
curve refers to the Chebyshev-Bezier curve over the interval [a, b] = [1,4] associated with the 
partition A, the red curve correspond to the partition (3, 2, 1), the green curve to the partition 
(4, 3, 2) and the blue curve to the partition (5, 4, 3). 



as a shape parameter. It would, therefore, be interesting to study the effect 
of standard operations on a fixed partition A, such as adding a box, removing 
a box, adding a row or cohimn and so on, on the shape of the curve. The 
problem is rather challenging and we will content ourself, here, with a simple 
experimental example. In Figure [3l we show the effect of adding boxes to the 
first row of the partition A = (2, 1). Adding boxes to the first row seems to have 
the effect of making the curve more and more far from the control polygon. 
However, adding the same number of boxes to every column seems to have the 
opposite effect as shown in Figure H) 

We can also define the tensor-product surfaces based on the Chebyshev- 
Bernstein basis associated with two different partitions. Namely, we can define 
a surface Tx,f_i by the parametric equation 






(81) 



where A and fi are partitions of length at most n, Pij are points in MP and 
(s, t) G [a, b] X [c, d]. Figure [5] shows an example of surfaces obtained from (|8T|). 
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Figure 5: Tensor-product surfaces obtained using equation II8II I: the transparent surface is 
associated with the parameters A = (2, 1), fj, = (1, 1) and [a, b] = [c, d] = [3, 4], while the blue 
surface is associated with the parameters A = (5, 1), fi = (5, 1) and [a, b] = [c, d] = [1, 6] 



.1. Derivative the Chehyshev- Bernstein Basis 

Consider a partition A = (Ai, A2, ..., A„) of length at most n in which we 



assume that 



Ai — A2. 



(82) 



Under the condition ([82|) . the derivative of the Chebyshev-Bernstein element 
BJ.^ is an element of the Chebyshev space E^m (n — 1) where A^"-* is the bottom 
partition of A. Therefore, the derivative of B^ ^ can be written as linear com- 
bination of the Chebyshev-Bernstein basis of Sxm {n — 1). Using the vanishing 
property of the Chebyshev-Bernstein bases stated in Theorem [TJ we derive as 
in (p5)) the following relationship 



dt 



un (k) 
^k.X 



(a) 



B, , ,,n. (a) 



^fc-l,A(n) (^) 



on (n-fc) 
^k,\ 



(b) 



fc-l,A(0) 



B 



_l (n-fc-l) 



fc.A(o) 



(6) 



Br.Ut)- (83) 



B" 



0. Let us denote 



in which we adopt the convention that i?"^ ^(o, = j^^^ ^(q, 

by 77 the partition rj = (A3, ..., A„). The partition 77 is the bottom partition of 

A'-^' and can well be written as ry = (A*^°^) , but for simplicity we will refer to 
this partition as 77. Inserting in (j83p the value of the derivatives from Theorem 
[SI we find 



B^iJ'Ha) 



n fx{n + l)^(n - 1) 5;,,o) (a"-^ &'=)^;,(o) (« 



n+l-k 



^fc-1) 



TDTl-l 

-°fc-l,A(0) 


(fc- 


(a) 


and 






S.",A^"- 


-k) 


(&) 



& - a /a(o) (")^ 



5'^(an+i-fc^ 5fe)5'^(a«-fc^ &fe-i) 



n h{n + l),^(n - 1) S^m (a"-^ 6'=)5^(o, (a 



n~-l~k 



5fe+i) 



R«-l (""'''"1)^(,\ 



b — a 



f\m (ny 



S'A(a"^^ &'=+i)5^(a"-fe-i, fofe) 



To write the formula for the derivative in a compact form, we define 

Definition 7. Let A = (Ai, A2, ..., A„) be a partition of length n > 2. Let a, b be 
real numbers. For k — 0, ...,n, we define the factor 

,0 ^ )bx(o-)(a ,b ) 



Rx{k,n) 



S\{o) (g" 

S'A(a"-^ &'=+i)S'„(a"-'=-i, 6'=) 



(84) 
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where A'"' — (A2,...,A„ is the bottom partition of X and ?/ = (A3,...,A„ is the 
bottom, partition of A*-"-' . 

From the last definition and by noticing that 

fxin + l)f^{n-l) _ 1 



we obtain 

Theorem 9. Let A = (Ai, A2, A3, ..., A„) be a partition of length at most n such 
that Ai = A2. Then the derivative of the Chebyshev-Bernstein basis associated 
with the partition A over an interval [a, b] satisfies 



A 



dt ib-a){ 



(i?A(fc - 1, ")B,ri.A(o, it) - Rxik, n)B'^-l„, it)) , 



A(o);„-i 

where X^*^' is the bottom partition of X and R\ik,n) is defined in equation l^84\ ) 
and in which rj is the partition rj = (A3,...,A„). We adopt here the convention 
that 5^7 ^^(0) = ^",A(o) = 0- 

Note that in the case the partition A is the empty partition, we recover 
the classical formulas for the derivative of the polynomial Bernstein basis. Let 
A = (Ai, A2, ..., A„) be a partition of length at most n such that Ai = A2 and 
Consider a fA("-)-function P written in the Chebyshev-Bernstein basis over an 
interval \a, b] as 



P(i)=^i?,"^,(t)P,. 



k=0 

Using Theorem |6j we can express the derivative of the function P in term of 
the Chebyshev-Bernstein basis over the interval [a, b] of the space E^m in — 1). 
Doing so leads to the following 

Theorem 10. Let X — (Ai, ..., A„) be a partition such that Ai — A2 and let P 
be a £\in)- function, written in the Chebyshev-Bernstein basis as 



P{t)=Y^Bl^it)P,. 

k=0 

Then, we have 

n-l 

where APi — -Pi+i — Pi and R\ik, n) is given in (f^l). 
In particular, we have 

na^^ A(n+l)g,(o)(a"-i,6) ^^ 

b-a fx(o)in) i)A(a",b) 

nb^^ /A(n + l)^;,(o)(Q,6"-i) ^^ 
^ i°) ~ 1 J r~\ c~r — T^\ — i^"^ ^ Pn-i) 



46 



(85) 



As from TheoremlHl we know the derivatives (i3";^)'(a) and (-B"_i a)'(^) ^^'^ ^y 
using the fact that the segments [Pq, Pi] and [P„_i, P„] are tangents to the curve 
at the point P{a) and P{b) respectively, we can show that the equations ([55)) 
are true independently if the partition A has it first two parts equals or not. 
It is rather interesting that computing the derivative {BQ^y{a) or {B'^^y{a) 
using the explicit expression of the Chebyshev-Bernstein basis pS)) reveal to be 
difficult. Equations (|55|) can be used to achieve the C^ continuity between two 
Chebyshev-Bezier curves associated with two different partitions, as follows 

Corollary 6. Let X = (Ai, A2, ..., A„ and /i — (/ii,/i2, ...,/x„) be two partitions 
of length at most n and let V = (Pq, Pi, ..., P„) (resp. (Q — Qo,Qi, ■■■,Qn)) 
be the Chebyshev-Bezier points of a £\{n) (resp. a £^{n)) function over an 
interval [a, b\ (resp. [&, c] j. Then the composite curve 7 formed by the two 
curves associated with the two control polygons V and Q is C^ at b if and only 
if Pn = Qo and 

nb^^ fx{n + l)S^^o,{a,b"-') nb^^ f^jn + 1) S^m{c,b"-') 

b-a f,,„,{n) S,{a,b-) ^^" ^-^^ " c - b /^,o, (n) S,.(c,6") ^^^ '^'> ' 

Example 7. Consider the Chebyshev-Bezier curve Fi of order n associated 
with the partition A = (1*^) with k < n and control polygon (Po,Pi, ...., Pn) 
over an interval [a,b]. Consider another Chebyshev-Bezier curve T2 of order n 
associated with the empty partition and control polygon (Qo, Qi, ■■■■Qn) over an 
interval [b, c] . From equations i85\) . a necessary and sufficient condition for the 
two curves Fi and F2 to be C^ at the point P„ is that 

Pn-Qo and ^^^^^^(i^.-^.-i)-^(Qi-Qo). 
k{b — a) ek(a,b"-) c — b 

If we denote by p the positive number such that P„ — P„_i = p{Qi — Qo), then, 
from the last equation, in order to achieve the C^ continuity at P„, we should 
choose the number c as 

u , fc(6-a)efc(a,b") , 

C=b-\--, 7— ; ; ; r- . (8d) 

(n+l)6pefe_i(a,&"-i) ^ ' 

Figure [^ shows the case n = 3 in this example, while Figure [^ shows another 
example of the application of Corollary \^ with the partitions A = (2,1) and 
M=(l,l). 



Remark 10. If a partition A = (Ai, A2, ..., A„) of length at most n satisfies 
Ai = A2 = ... = As, s < n, then we can iterate Theorem [S] to compute the 
derivatives up to order s — 1 of the Chebyshev-Bernstein basis. 

Consider, now, a partition A = (Ai, A2, ..., A„) of length at most n in which 
we assume this time that Ai 7^ A2. Under this condition, the derivative of 
the Chebyshev-Bernstein element B^ ^ over an interval [a, 6] is an element of 
the Chebyshev space £^(n) where /i is the partition p = (Ai — 1, A2, ..., A„). 
The derivative of B^ ^ can be written as linear combination of the Chebyshev- 
Bernstein basis of £fj,{n). However, if we use the vanishing properties of the 
Chebyshev-Bernstein bases, we arrive to a three-term recurrence relation be- 
tween the two Chebyshev-Bernstein bases and in which Theorem |5] does not 
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Figure 6: C^ continuity at the point P3 between two Chebyshev-Bezier curves associ- 
ated with two differents partitions. The Chebyshev-Bezier curve with Chebyshev-Bezier 
points (Pq^ Pi, P2^ P'i) is associated to the partition (1, 1) and parametrized over the interval 
[a,b] = [1,3]. The Chebyshev-Bezier curve with Chebyshev-Bezier points (Qoi Qli Q2, Q3) is 
associated to the empty partition and parametrized over the interval [3, c], where the param- 
eter c was computed using equation II86I I to achieve the C^ continuity, (see example 13) 




Figure 7: C'^ continuity at the point P3 between two Chebyshev-Bezier curves associ- 
ated with two differents partitions. The Chebyshev-Bezier curve with Chebyshev-Bezier 
points {Po, Pi, P2, P'i) is associated to the partition (2, 1) and parametrized over the interval 
[a, b] = [1,3]. The Chebyshev-Bezier curve with Chebyshev-Bezier points (Qo, Ql, Q2, Q3) is 
associated to the partition (1, 1) and parametrized over the interval [3,c], where the parameter 
c was computed using the conditions of corollary [6] to achieve the C^ continuity. 
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allow for an easy way to compute the necessary coefficients. To solve this prob- 
lem, we can instead proceed as follows : As Ai 7^ A2, we have necessarily Ai > A2. 
Therefore, the partition rj = (Ai — 1, Ai — 1, A2, ..., A„) is a dimension elevation 
partition of A. We can compute the Chebyshev-Bernstein basis of £^(71 + 1) as 
a function of the Chebyshev-Bernstein basis of £\ (n) according to Theorem [T] 
Now the partition rj satisfy the condition that its first two parts are equals, and 
therefore, we can use Theorem [5] to compute the derivative. Proceeding along 
these two steps, in which we omit the computation as they can be readily done, 
we find 

Theorem 1 1 . Let A = (Ai , A2 , A3 , . . . , A„ ) be a partition of length at most n such 
that Ai 7^ A2. Then the derivative of the Chebyshev-Bernstein basis associated 
with the partition A over an interval [a, b] satisfies 

",, ' = ,, ^ ,7,^ (Gi(fc, n)Bl_,^,{t) + G2(fc,n)B^.„(t) + G3(fc,n)B,V,^(t)) , 
at (^t, - a){^) ^^ 

where 

Gi{k,n) = (n + I - k)aT\(n,k)R,^{k - l,n + 1), 

G2{k,n) = Rri{k,n + l)((fc + l)6A^(n, fc) - (n + 1 - fc)ar^(n, fc)), 

and 

G3(fc, n) = -{k + l)bAl{n, k)R^{k + 1, n + 1). 

and the partition rj and fi are given by rj = (Ai — l,Ai — 1,A2,...,A„) and 
ji — (Ai — 1, A2, ..., A„), the factors A^, T\{n,k) and i?^ are defined in \6T^ , 
\68\) and \84-^ respectively. We adopt the convention that ^"i_^ = B"+i,^ = 0. 

Let A — (Ai, A2, ..., A„) be a partition of length at most n such that Ai ^ A2 
and consider a £"> (n)-function P written in the Chebyshev-Bernstein basis over 
an interval [a, b] as 

n 
fe=0 

Using Theorem llll we can express the derivative of the function P in term of 
the Chebyshev-Bernstein basis over the interval [a, b] of the space £^{n). Doing 
so leads to the following 

Theorem 12. Let A — (Ai, ..., A„) be a partition such that Ai ^ A2 and let P 

be a E\{n)- function, written in the Chebyshev-Bernstein basis over an interval 

\a, 61 as 



P(i)^^i?If,(t)P,. 



fc=0 



Then, we have 



v\ 



where the factors Gi, G2, G3 and the partitions rj and /i are defined in Theorem 
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9. Conclusions 

In this paper, we carried out a comprehensive study of the notion of Cheby- 
shev blossom in Miintz spaces, thereby showing their adequacy in free-form 
design schemes. An interesting aspect of the work was the followed method- 
ology in providing for an explicit expression of the Chebyshev-Bernstein basis. 
Most of the steps in the proof were combinatorial in nature. Similar arguments, 
therefore, could be applied to any extended Chebyshev space constructed from 
weight functions, in the sense that the condensation formula will provide us 
with the pscudo- affinity factor and in which the combinatorics of the de Castel- 
jau paths can be employed to give extra-information on the derivatives of the 
Chebyshev-Bernstein basis. Such a program will be the object of a forthcoming 
contribution. Moreover, the problem of higher order continuity and the issue of 
shaping with Young diagrams lead to interesting problems for future work. 

Acknowledgement : This work was partially supported by the MEXT Global COE 
project. 
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